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I . INTRODUCTION 


i 

I 

I 

I The  Air  Transport  of  Radiation  (ATR)  code  has  been  de- 

veloped to  provide  a fast  running  computer  code  to  predict 
t energy- angular  radiation  environments  from  neutron,  gamma-ray 

and  X-ray  sources  in  the  atmosphere.  During  the  development  of 
ATR,  there  have  been  several  reports describing  the  capa- 
bilities of  each  version  of  the  code  and  the  generation  of  the 
data  bases  on  which  the  parametrizations  were  based.  This  re- 
port summarizes  the  total  capability  of  ATR-4  - perhaps  the  most 
definitive  time- independent  version  of  ATR.  Detailed  descrip- 
tions of  the  generation  and  parametrization  of  each  data  base 
component  were  described  in  the  previous  reports  and  only  a brief 

description  is  provided  here.  Emphasis  is  placed,  however,  on 
the  capabilities,  limitations  and  use  of  the  code.  As  such,  this 
report  may  be  considered  as  a user's  manual  for  the  ATR-4  code 
' and  except  for  special  cases  where  the  ATR  code  has  been  cus- 

tomized to  a particular  computer  capability,  this  version  of 
the  code  should  replace  all  earlier  versions.  The  time-depen- 
dent version  TDATR  is  documented  elsewhere  and  is  treated  as  a 
separate  code. 

The  Defense  Nuclear  Agency  (DNA)  has  for  some  time 
sponsored  efforts  to  develop  sophisticated  transport  techniques 
and  to  accurately  evaluate  cross  sections  in  order  to  increase 
the  confidence  in  and  precision  of  atmospheric  radiation  trans- 
port calculations.  ATR  was  designed  to  embody  the  results  of 
I these  efforts  and  make  them  available  to  the  general  community 

I in  a form  that  eliminates  the  need  for  data  tables  and  hand 

1 

; folding.  ATR  does  not  perform  radiation  transport  calculations 

I but  contains,  instead,  parametric  representations  of  fluence 

I results  derived  from  such  calculations.  This  approach  enables 
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the  ATR  users  to  obtain  satisfactory  results  without  requiring 
either  radiation  transport  expertise  or  programming  sophistication. 
Figure  1 illustrates  the  characteristics  of  the  ATR  model. 


Figure  1.  ATR  Model 
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A summary  of  the  development  of  the  ATR  code  is  presented 
in  Table  1 which  describes  reports  which  have  been  published 
with  a very  brief  description  of  their  content.  Since  ATR  has 
been  an  evolving  code,  changes  in  later  versions  have  made 
the  previous  ones  incomplete  or  obsolete.  In  addition,  new 
features  such  as  a simple  fission  product  model  have  been  added 
and  are  documented  in  this  report.  Thus,  this  report  serves  to 
document  the  new  features  in  ATR-4  as  well  as  to  summarize  the 
capabilities  of  the  complete  code  package. 


TABLE  1. 


ATR  DOCUMENTATION  SUMMARY 


Report 


Content 


SAI-71-565-LJ 
(November  1971) 


DNA  28031 
(May  1972) 


DNA  314AZ 
(April  1973) 

DNA  3279T 
(August  1974) 

DNA  3362Z 
(August  1974) 

DHA-38I9F 
(July  1975) 


DNA  4061 
(January  1976) 


Paper  given  at  RSIC  Workshop  on  Radiation 
Transport  in  Air.  Describes  logic  of  code 
and  some  data  ba.se  development. 

First  report  describing  basic  concepts  of 
ATR  and  the  first  distributed  version  of 
the  code- - describes  neutron  and  secondary 
gamma  ray  data  base  generation.  Includes 
air-ground  and  exponential  'ir  correction 
factors . 

Users  manual  for  ATR-2  version  of  the  code-- 
does  not  describe  data  base  generation. 

Describes  data  base  generation  for  photons 
(prompt  gamma  rays  and  X-rays)  and  the  prompt 
gamma  ray  air-ground  correction  factors. 

Summary  of  the  capabilities  of  the  ATR  code 
with  updates  to  ATR-2. 

Describes  ATR-3  including  new  data  base  using 
DNA  cross  section  library,  new  air-ground  cor- 
rection factors,  low  energy  X-rays  and  new 
REGROUP  routine. 

Describes  TDATR,  the  time  dependent  prompt 
photon  and  secondary  gamma  ray  version  of  ATR. 


This  report 


Describes  fission  product  model  and  summarizes 
total  capability  of  the  ATR-4  code. 


2.  GENERAL  DESCRIPTION  OF  THE  ATR  CODE  SYSTEM 


The  ATR  code  was  developed  to  provide  an  easy-to-use 
means  of  determining  the  complete  radiation  environment  from 
a nuclear  weapon  burst  in  the  atmosphere.  An  "engineering"  ap- 
proach was  taken  in  that  the  results  of  detailed  calculations 
were  parametrized  to  provide  a compact  data  base  instead  of 
providing  a table  of  results  from  a radiation  transport  code. 
This  "engineering"  approach  results  in  a trade-off  in  accuracy 
and  data  storage  requirements.  In  addition  the  treatment  of 
the  effect  of  the  air-ground  interface  and  exponential  atmo- 
sphere is  provided  by  utilizing  "correction  factors" 
rather  than  parametrizing  many  detailed  transport  results. 
Because  of  the  lack  of  detailed  data  bases  for  these  two  ef- 
fects, the  "engineering"  approach  may  be  inadequate  for  some 
applications,  but  the  philosophy  of  utilizing  the  best  data 
available  in  a consistent  manner  has  been  adopted.  Due  to 
lack  of  data  available,  correction  factors  are  only  applied 
to  integral  quantities. 

The  ATR  code  is  a user  oriented  code  which  can  be  executed 
in  either  a batch  mode  or  from  a desk  teletype  unit.  For  many 
problems  the  built-in  sources  and  responses  may  be  adequate, 
thus  reducing  the  amount  of  input  data  that  is  required.  Yet 
the  flexibility  to  investigate  a very  specific  case  is  main- 
tained by  allowing  all  aspects  of  the  problem  to  be  input. 

Table  2 lists  the  components  of  radiation  that  can  be  calculated 
by  ATR.  Results  are  provided  at  fixed  detector  locations  with 
any  source,  detector  or  ground  elevation  allowed.  Figure  2 
shows  the  definition  of  the  geometry  models  utilized.  The 
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TABLE  2.  RADIATION  ENVIRONMENTS  CALCULATED  BY  ATR 


Types  of  Radiation: 

Neutrons 

Secondary  Ganrnia  Rays 
Prompt  Gamma  Rays 
X-rays 

Fission  Product  Gamma  Rays 


Results : 


Energy- angular  fluence  and  current 
(except  for  fission  product  gamma  rays) 

Angle  integrated  results 

Energy  weighted  results 


mass  thickness  of  air  is  determined  between  the  source  and 
detector  and  results  are  scaled  to  an  equivalent  distance  in 
homogeneous  air.  The  data  base  was  calculated  with  an  air  den- 

3 

sity  of  1.11  mg/cm  . 

The  fundamental  calculation  in  ATR  is  to  determine  the 
free-field  fluence  for  a given  problem  geometry  and  source 
spectrum.  In  order  to  calculate  such  problems  for  the  actual 
atmosphere,  two  assumptions  are  made: 

1.  The  transport  of  radiation  between  source  and 
target  can  be  approximated  by  homogeneous  air 
transport  results  for  a density  equal  to  the 
average  density  between  the  source  and  target 
altitudes . 

2.  The  corrections  due  to  air-ground  interface 
and  exponential  air  density  effects  can  be 
included  as  perturbations  to  homogeneous  air 
results . 
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~1035  g/cm^  VACUUM 


Non-Uniform  Air  Corrections 
~913  g/cm^  15,000  Meters 
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Figure  2.  ATR  Geometry  Model 
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The  general  calculational  scheme  in  ATR  is  therefore  to  scale 
transport  results  from  the  parametric  model  density  of  1.11 

3 

mg/cm  to  the  average  density  of  the  problem  and  then  to  in- 
clude correction  factors  for  the  effects  of  an  air-ground 
interface  and  exponential  air  density.  These  correction  factors 
were  essentially  negligible  for  altitudes  from  1-15  km.  Density 
scaling  laws  for  the  Boltzmann  equations  show  that  the  angular 
flux  ^)^(E,J^,r)  in  the  transformed  system  of  density 
related  to  the  flux  in  the  original  system  at  density  p(r)  by 

(|>(E,n.r)  = !j>^(E,fi,r^) 


where 


and 


i . e . , for 


r^  = Kr 


Pj.(r^)  = p(r)/K 


Pt’^t  = 


This  is  the  well  known  result  that  the  flux  depends  only 

on  the  mass  of  air  through  which  the  neutrons  and  photons  are 

2 

being  transported.  The  ATR  parametric  model  is  based  on  U-nr 

2 

quantities  as  a function  of  penetration  distance  in  gm/cm 

The  total  mass  of  air  between  source  and  target  for  a given 

_ 2 

problem  is  calculated  from  R p and  used  to  obtain  4ttR  quan- 

s s 

tities  from  the  parametric  data  base,  i.e., 


AttR" 


) = AttR  (J) 

s o o 


For 


= R p 
o o 
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where 


p = average  air  density  of  the  problem 
pQ  = air  density  of  ATR  data  base  (1.11  mg/cm^) 

~ radiation  free-field  fluence  calculated  by  para- 
metric  equations  at  density  p^  and  distance 

4)  = radiation  free-field  fluence  of  the  geometry  of 

® interest. 

The  average  air  density  of  the  problem  is  calculated  from 


where  p(Z)  is  the  variation  in  the  atmospheric  density 
as  a function  of  altitude  above  sea  level. 

The  atmospheric-density  is  modeled  from  the  1962  U.S. 
Standard  Atmosphere  for  altitudes  from  0-300  km  based  on  a 
product  of  32  exponential  terms  and  is  used  for  altitudes  from 
0-100  km.  Above  100  km,  an  empirical  curve  fit  is  used. 

The  data  base  of  ATR-4  was  calculated  for  neutrons,  second- 
ary gamma  rays,  prompt  gamma  rays  and  X-rays  using  the  ANISN^^^ 
discrete  ordinates  code  and  utilizing  the  DNA  cross  section  li- 
brary^  ' for  all  but  X-rays.'  ' Results  as  a function  of  source 
energy  and  range,  detector  energy  and  angle  for  ranges  to  550 
gm/cm  were  deteinnined.  Since  the  earlier  versions  of  ATR  were 
based  on  neutron  and  secondary  gamma-ray  results  determined  with 
the  cross  sections  evaluated  by  Straker,  the  new  data  base  was 
produced  by  the  parametrization  of  correction  factors  to  the 
earlier  data  base.  Thus  it  is  possible  to  evaluate  differences 
due  to  cross  sections  since  the  earlier  data  base  was  not  de- 
stroyed. The  fission  product  gamma-ray  capability  in  ATR-4  is 
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based  on  simple  empirical  models,  and  the  basis  of  the  radiation 
transport  used  in  the  fission  product  model  comes  from  pervious 
ATR  techniques.  A more  detailed  discussion  of  the  data  bases 
is  given  in  Section  3. 

One  of  the  main  features  of  the  ATR-4  code  as  well  as  pre- 
vious ATR  models  is  the  use  of  control  commands  for  selecting  and 
executing  options.  The  list  of  control  commands  for  setting  up 
an  ATR  run  is  discussed  in  Section  6.  In  general,  the  order  of 
control  commands  is  immaterial  for  any  given  set  of  problems . The 
initialization  routine  contains  some  calculations  on  parameters 
which  can  be  evaluated  on  the  basis  of  other  parameters  as  well 
as  some  storage  initialization.  When  a command  is  read  in  and 
interpreted,  the  corresponding  subroutine  is  called  to  interpret 
and  properly  place  all  of  the  parameters  corresponding  to  the 
command.  Errors  in  control  command  format  or  structure  are  noted 
immediately  but  no  other  visible  action  is  taken  by  the  code  until 
the  "EXC  command  is  used,  at  which  point  output  is  generated  which 
may  be  just  a few  numbers  or  hundreds  of  pages,  depending  on  the 
input  parameters. 

The  code  also  contains  a set  of  subroutines  which  permit 
an  essentially  format-free  placement  of  numbers.  The  reading 
of  values  is  upward  compatible  in  that  real  values  need  not 
contain  decimal  points.  Control  commands  and  their  associated 
parameter  values  are  cumulative  so  that  even  after  the  results 
of  a particular  problem  have  been  displayed,  the  parameter 
values  and  flags  remain  the  same  for  the  next  problem  unless 
specified  anew  with  another  command  of  the  same  type.  Con- 
sequently, only  the  most  essential  parameters  need  to  be 
specified  from  one  problem  to  the  next.  This  feature  allows 
multiple  runs  to  be  made  easily  which  require  only  a few  para- 
meter changes  to  be  made  for  each  problem.  In  order  to  clear 
out  unwanted  parameters  the  *STOP  command  is  used  after  which 
the  problem  must  be  completely  respecified. 
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2.1  r i'i.METRIZATION  ACCURACY  CONSIDERATIONS 


;g  the  accuracy  of  ATR  parametrization  of  one 
dimens  ^tion  transport,  the  convention  is  adopted 

that  the  dw.Ui.av.y  is  relative  to  the  data  base  from  which  the 
fits  were  produced.  Comparisons  with  other  independent  calcu- 
lations using,  say,  different  cross  sections,  energy  group 
structure,  or  transport  methods  should  first  be  made  with  the 
data  base.  In  this  manner,  the  accuracy  of  the  parametriza- 
tion can  be  separated  from  the  problems  associated  with  making 
comparisons  of  independent  transport  calculations.  The  data 
base  used  for  parametrization  of  the  neutron  and  secondary 
gamma  ray  transport  is  similar  to  the  data  reported  by  Straker 
and  Gritzner. The  generation  of  the  data  bases  for  the 
prompt  gamma  rays  and  X-rays  are  described  in  a document  by 
Woolson,  et  al . ' 

The  parametrization  of  the  data  bases  contains  three 
separate  fluence  quantities  which  were  independently  fit.  These 
quantities,  the  total  fluence,  the  normalized  scalar  fluence, 
and  the  angular  dependent  ratio  function  are  described  in  fol- 
lowing sections.  During  the  fitting  process,  accuracy  standards 
were  established  for  each  of  the  fluence  quantities.  These 
standards  are  given  in  Table  3.  The  parametrization  was  based 
on  finding  techniques  which  would  fit  most  of  the  data  to  the 
average  accuracy  standards  given  in  the  table.  A second  order 
correction  method  was  then  employed  to  improve  the  accuracy  of 
the  fits  which  exhibited  anomalous  errors. 

The  final  parametrizations  incorporated  into  the  ATR  code 
are  generally  accurate  to  within  the  stated  standards.  In  parti- 
cular, the  total  fluence  functions  and  normalized  scalar  fluence 
were  found  to  be  generally  much  more  accurate  than  the  standards 
indicate.  The  ratio  function  fits,  which  contains  the  informa- 
tion about  the  angular  variation  of  the  fluence,  were 
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TABLE  3.  ACCURACY  STANDARDS  FOR  PARAMETERIZATIONS  OF 
INFINITE  HOMOGENEOUS  AIR  TRANSPORT  IN  ATR 


Fluence  Type 

Definition 

Accuracy 

Standard 

Total  fluence 

T(Eg,R)  = //  (j)(Eg,R,E^,y)  dp  dE^ 

57o 

Normalized 

scalar 

fluence 

/ (t)(E  ,R,E,,y)  dp 


107o 

Ratio  function 

4i(E  ,R,E  ,,p) 

R(Es.R.Ed.y)  T(Eg,R)  F(Eg.R,E^) 

207o 

compared  only  on  the  average  with  the  data  base  values  for  prompt 
gamma  rays  and  X-rays  since  the  angular  fluence  calculated  by 
discrete  ordinates  for  photon  transport  usually  exhibit  non- 
physical variations  due  to  the  strongly  peaked  angular  distri- 
butions. Although  the  scalar  fluences  (integrated  over  angle) 
are  accurate,  the  differential  angular  fluence  required  some 
smoothing.  Thus,  the  comparisons  of  the  ratio  function  fits 
were  made  with  the  smoothed  angular  fluences  and  not  directly 
with  the  individual  data  base  values. 
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3.  ATR  DATA  BASES  AND  THEIR  PARAMETRIZATION 

Different  components  of  the  ATR  data  base  have  been 
developed  over  a period  of  about  five  years.  The  first  neu- 
tron and  secondary  gamma  ray  data  base  utilized  the  cross 
sections  assembled  by  Straker  around  1967.  This  data  base 
was  used  in  Versions  1 and  2 of  ATR  and  was  replaced  by  a new 
data  base  generated  with  DNA  cross  sections  evaluated  around 
1974.  The  new  data  base  for  neutrons  and  secondary  gamma  rays 
was  first  distributed  in  Version  3 of  ATR  and  because  of  the 
need  for  comparisons  with  older  calculations,  the  new  data  base 
was  incorporated  as  a modification  to  the  original  data  base. 

The  prompt  photon  data  base  first  appeared  in  Version  2 
of  ATR  with  modifications  to  the  X-ray  data  base  to  extend  it 
to  lower  energies  in  Version  3.  A delayed  gamma  ray  (fission 
product)  data  base  was  added  in  Version  4.  A separate  version 
of  ATR,  TDATR,  contains  the  time  dependent  prompt  gamma  ray 
and  secondary  gamma  ray  data  base. 

The  ATR  data  base  consists  of  four  separate  entities:  a 
neutron  and  secondary  gamma  ray  data  base,  a prompt  gamma  ray 
data  base,  an  X-ray  data  base  and  a fission  product  data  base. 

The  first  three  sets  of  data  were  generated  in  infinite,  homo- 

3 

geneous  air  at  a density  of  1.11  mg/cm  by  the  discrete  ordi- 
nates radiation  transport  technique.  The  transport  calculations 

2 

were  performed  out  to  a range  of  550  gn/cm  for  neutrons,  secondary 
gamma  rays  and  prompt  gamma  rays.  For  the  X-ray  calculations  the 
range  was  dependent  on  source  energy  but  covered  at  least  15 
mean-free-paths  for  uncollided  (source)  particles. 
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Each  radiation  transport  data  base  consisted  of  energy  and 
angular  dependent  fluence  as  function  of  range  for  individual 
source  energy  bands.  Thus,  each  separate  data  base  can  be 
described  in  terms  of  a four-dimensional  array  of  fluence 
values,  defined  by 


-f 


k+l 


’'"j+l 


^sikj  " J f du  |)  g(r^  ,E,  u) /AEj^ 

*^k 


where 


p “ four-dimensional  array  of  fluence  values  com- 

® ^ puted  by  discrete  ordinates  with  source  energy 

index,  s;  range  index,  i;  receiver  energy  index, 
k;  and  angle  index,  j;  for  neutrons,  secondary 
gamma  rays,  prompt  gamma  rays,  and  X-rays.  The 
units  on  the  fluence  quantity  are  particles 
(neutron , photon)  /cm2-MeV''-steradian-source  "par- 
ticle" in  energy  group  s. 

^ ,(r.,E,u)  = radiation  fluence  environment  at  range  r^ , 

^ ^ energy  E and  angle  cosine  p in  particles/ 

cm2-MeV-steradian  per  unit  source  in  energy 
group  s . 

r.  = range  mesh  point,  in  air  at  density  1.11 
^ mg/ cm3. 

(Ej^, = lower  and  upper  limits  of  energy  band  k. 

" \+l  ■ \- 

= lower  and  upper  limits  of  angle  cosine  bin  j. 

Awj  ■ - Wj. 

In  the  following  the  average  value  of  the  energy  within  a 
group  denoted  by  E^  and  defined  by 


keV  for  X-rays. 
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J 


r 


T 


E 

n 


E +E  , , 
n n+1 


is  used.  Similar  relationships  will  be  used  for  the  angular 
bins.  Additionally,  the  notation  Ej^,  E^^  will  be  used  to 
indicate  the  lower  and  upper  energy  limits  of  a particular 
data  base. 

The  ATR  m.odels  of  the  radiation  environment  retain  the 
group-wise  features  of  the  data  base  for  source  energy  groups, 
receiver  energy  groups,  and  angle  cosine  groups  but  reconstruct 
the  environment  as  a continuous  function  of  r.  The  ATR  scattered 
fluence  model  was  constructed  as  a product  of  three  functionals 
which  exhibited  slower  variation  as  a function  of  range  than  the 
basic  fluence  quantity.  The  fluence  model  including  the  un- 
collided contribution  is  given  by 


-or 

s 


'"sk 

• "Set 

k 


6 . . 

JJi 


27iAy . 

Jh 


where 


Ts(r) 


= ATR  fluence  model  at  range,  r,  for  source  band 
s;  receiver  band,  k;  and  angle  cosine  band,  j. 
Units  are  particles/cm^-MeV-steradian-source  in 
group  s . 

= total  scattered  fluence  in  units  of  particles/ 


cm 


2 _ 


source  in  group  s. 


''skj  <'•> 


c 


normalized  scalar  fluence  shape  in  units  of 
{fraction  of  Tg(r)  in  energy  group  k)/MeV. 

ratio  function  in  units  of  (fraction  of 
Tg(r)  • Fsk(r)  in  cosine  bin  j } /steradian. 

total  macroscopic  multigroup  cross  section  for 
energy  group  s. 


ik 


= kronecker  delta  function. 


j,  = forwardmost  angular  bin, 


k.. 
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In  this  representation,  most  of  the  spatial  variation  of 

the  fluence  is  built  into  the  total  fluence  function  T (r) ; 

s 

the  energy  dependence  is  contained  in  the  scalar  fluence  shape 
sk' 


F_,_(r);  and  the  angular  distribution,  in  the  ratio  function. 


Several  relationships  exist  between  the  real  environment 
4)g(r,E,M),  and  the  ATR  modeled  fluence  • These  relation- 

ships are  summarized  in  Table  4. 


TABLE  4.  SCATTERED  FLUENCE  RELATIONSHIPS 


r 

ATR  Model 

Radiation  Environment 

1 

! 

^skifrl  ; 

1 

i 

®^k+l  “j+l 

y dE  y da0g  (r,  E, /i)  Aej^  A/J. 

E,  u . ' 

k 3 

Tg  fr)  2n  J dE  J (r,  E, 

1 E,  -1 

I ^ 

^sk 

2n  J AE  J da0g  (r,  E,  ^i) 

2,rAE^  f J (r.  E,  a) 

«skj 

y (r,  E,  ft) 

K 1 1 

k+1 

2nAft^  y dE  f (*■>  E,  ^‘) 

\ -1 

Notation  defined  in  the  text,  uncollided  fluence  not  included. 
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Note  that  with  the  fluence  model  utilized  by  ATR  there 
may  be  several  ways  to  construct  a given  fluence  quantity.  For 
example,  the  energy- dependent  scalar  fluence  could  be  constructed 
from  the  ATR  models  by  the  operation 


however,  it  can  be  also  retrieved  by 


-o  r 
s 


T (r)  F , . . + 
s^  sk(r) 


Airr^ 


The  second  alternative  is  more  accurate  because  it  does  not 
contain  possible  errors  present  in  the  ratio  function  fit.  Inter- 
nally ATR  employs  the  most  accurate  representation  in  calculating 
fluence  quantities. 

The  ATR  code  always  performs  integral  mass  scaling  to 
correct  for  density  differences  between  the  input  problem 
geometry  and  the  data  base  calculation.  If  necessary,  cor- 
rections are  made  for  air-ground  interface  or  high  altitude 
perturbations.  The  resulting  scaled  and  corrected  fluences, 
denoted  by  the  primed  quantities,  are  functionally  related  to 
infinite  homogeneous  models  by 


skj 


(r)  = 


where  h^  is  the  ground  height  relative  to  sea  level; 
are  the  source  and  target  heights  respectively  relative  to 
the  ground  height;  and 


y 


h +h. 
' g t 


1. 11x10’^  l^s'^t ' K'„+h 

g ° 


p(h)  dh 
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where  p(h)  is  the  atmospheric  density  at  height  h above 
sea  level.  Detailed  information  concerning  the  scaling  con- 
siderations can  be  found  in  Ref.  1.  The  ATR  output  quan- 
tities will  now  be  given  in  terms  of  operations  on  the  scaled 
and  corrected  data. 

Energy-Angular  Dependent  Fluence 


particles 


s — T 

cm  "MeV • steradian 


where 


S 

s 


source  strength  in 


source 

MeV 


for  energy  band  AE^. 


Angular  Dependent  Fluence 


Fjr)  = I I 


^ s k ^ “ cm^  • steradian 


particles 


where  AEj^  is  the  receiver  energy  band  for  group  k. 
Energy  Dependent  Fluence 


Fu(r)  = I S. 


-o  r 
s 


T;(r)  F;^(r)  AE^  + 


Arr 


T-  %k 


particles 

2 

cm  -MeV 


Total  Fluence 


F(r)  - I S 
s 


T' (r)  + 


-o  r 
s 


particles 

—2 

cm 


Energy-Angular  Dependent  Energy  Fluence 


MeV 


s ^ skj 


cm  ’MeV* steradian 
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Angular  Dependent  Energy  Fluence 


e.  (r)  = I I (r)  E,,  AE,,  AE 


MeV 


s k 


k k s 2 

cm  'steradian 


Energy  Dependent  Energy  Fluence 


- I S,  E 

s 

Total  Energy  Fluence 


-o  r 

H-  «sk 

4TTr 


MeV 

T 


cm  'MeV 


e(r)  = 5;  ^ S E,  AE 
s k ® ® 


-a  r 

r (r)  F;^(r)  AE^  + ^ 


MeV 

~T 

cm 


Energy-Angular  Dependent  Number  Current 

"kj  ■ ps  ♦ ^kj  "j  “s  -7 

s cm  ‘MeV* steradian 


where  and  the  cosine  boundaries  of  cosine  group 


j,  and  y^  = (yj+yj_^^) /2 . 


Angular  Dependent  Number  Current 


cm  'steradian 


Energy  Dependent  Number  Current 


Ak<')  = U S,  2"  AW,  AE^  Earticles 

S J J J J s 


Total  Number  Current 


t(r)  = I I I 2TT  Ay,  AE,  AE„ 


s j k 


J k s 


cm 
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Energy-Angular  Dependent  Energy  Current 


e>,  . (r)  = y S . (r)  E,  y.  AE 

kj  ^ s skj  k j s 


cm  'MeV* steradian 


Angular  Dependent  Energy  Current 


cm  • steradian 


Energy  Dependent  Energy  Current 


c,.(r)  = y I S E,  y.  2tt  Ay  AE  — ^ 

k . s stcj  k j J s 


Total  Energy  Current 


c(r)  = II  I \ 


Energy-Angular  Dependent  Response 


- J =s  \ HeV-lgerldlan 


where  Rj^  is  the  response  function  in  units  of  response/ 
particles/cm2  for  energy  group  k. 

Angular  Dependent  Response 


i '»;kj<’^)  \ “s  stiradlan 


Energy  Dependent  Response 


-o  r 
s 


M - ^ Sa  \ '';k<’^>  ^ «sk 


response 


W 


Total  Response 


R(r)  = I I R, 
s k 


AE 


-o  r 
s 


T;(r)  F;j^(r)  + 


Anr 


sk 


response 


2 

An  option  exists  in  ATR  to  obtain  "Anr  " quantities.  When 

2 

Airr  " output  is  desired,  the  previous  defined  quantities  are 

2 2 

multiplied  by  the  Anr  constant  and  units  of  cm  are  placed 
in  the  numerator  of  the  unit  definition. 


3 . 1 NEUTRONS  AND  SECONDARY  GAMMA  RAYS 

The  data  base  for  neutrons  and  secondary  gamma  rays  was 
generated  using  the  DNA  cross  sections  for  nitrogen  and  oxygen 
as  they  existed  in  early  197A. 

The  cross  section  data  actually  used  in  generating  the  ATR 
data  base  were  in  the  multigroup  form.  Briefly,  multigroup  cross 
Sections  are  produced  by  expressing  the  angular  variation  of  the 
scattering  cross  section  as  constants  in  a Legendre  series  ex- 
pansion and  performing  weighted  averages  of  all  cross  sections 
over  the  energy  variable.  The  AMPX^^^^  code,  which  at  the  time  was 
only  fully  operational  at  the  Oak  Ridge  National  Laboratory  (ONRL) , 
was  used  at  ORNL  to  generate  the  multigroup  cross  sections.  The 
same  22  neutron  group  and  18  photon  group  energy  structure  in 
previous  ATR  versions  was  used.  A hand  modification  was  made 
to  add  10.8  MeV  nitrogen  capture  gamma-rays  to  the  top  photon 
group  (8-10  MeV).  The  cross  sections  used  in  generating  the 
new  ATR  neutron  and  secondary  gamma-ray  transport  data  base  can  be 
obtained  from  ORNL.  Radiation  Shielding  Information  Center  (RSIC) . 

The  radiation  transport  calculations  were  performed  us- 
ing the  ANISN  one-dimensional  discrete  ordinates  code.  The 
calculations  were  carried  out  in  a manner  similar  to  the  pre- 
vious data  base  calculations.  Gaussian  quadrature  of  order 
16  and  a Legendre  expansion  were  employed.  The 
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calculations  were  made  in  spherical  geometry  out  to  a range  of 
5 km  in  air,  of  density  1.11  mg/cm'  and  185  radial  intervals 
were  used.  A separate  transport  calculation  was  performed  for 
a uniform  source  in  each  of  the  top  18  neutron  energy  groups. 


If  ifi  g (Eg  , , r , m)  represents  the  old  data  base  for  source 
energy  group  E^ , detector  energy  group  E^,  range  r and  angle 
u and  ‘J)  ^(Eg  , E^ , r , vO  represents  the  new  data  base  with  the 
latest  DNA  cross  sections  then  by  integrating  both  quantities 
over  angle,  the  following  ratio  function  is  obtained: 


R(E^.Ed.r) 


4’n(f'„>Ej,r) 


Previous  versions  of  ATR  made  available  the  approximation 
ii^^g(E  ,Ej,r)  to  the  Straker  data  base.  Thus  by  parametrizing 
R(E^!Ed.r) 

- *«(Ed.Ed.r)  *5(E^,Ed.r) 


is  the  representation  of  the  new  data.  In  actual  practice,  the 
following  ratio  function  was  parametrized: 


R(Eg,E^,r) 


^^g(Eg,E^,r) 


The  ratio  function  for  neutrons  was  relatively  straight- 
forward to  curve  fit  by: 

2 

R(E  ,E,  ,r)  = exp(a„  + a,r  + a„r  ) 

Sq  Oq  u j.  z 

that  is,  for  a given  source  and  detector  energy  group  the  log 

of  the  ratio  function  is  presented  as  a polynomial  function  of 

2 

distance  in  gm/cm  . Tiius , t he  parametrization  coefficients  are 
functions  of  source  and  detector  energy.  For  some  sources  even 
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the  quadratic  term  was  not  necessary.  Fewer  than  600  coefficients 
were  necessary  to  completely  represent  the  neutron  ratio  function. 

For  the  secondary  gamma  rays, the  curve  fitting  procedure 
was  significantly  more  complicated.  Each  curve  was  examined 
as  a function  of  distance  and  the  more  complicated  curves  re- 
quired fourth  order  polynomials.  Furthermore,  there  were  some 
ratio  functions  which  required  a function  of  the  following  form: 

R(Eg,E^,r)  = exp(aQ  + a^r  + + a^r^  ^ ^5 

for  adequate  representation.  In  all,  nearly  1200  coefficients 
were  necessary  to  represent  the  secondary  gamma  ray  ratio  func- 
tions . 

A final  comparison  was  made  with  the  data  base  to 
make  sure  that  the  approximations  did  not  compound  in  an  un- 
reasonable manner.  The  spectral  comparisons  show  that  the 

maximum  error  of  4)r>(E  ,E,,r)  versus  .E,,r)  is  on  the 

usd  Usd 

order  of  10  percent. 

After  the  spectral  parametrization  was  complete,  com- 
parisons of  the  angular  dependence  were  made  and  large  errors 
were  "massaged"  or  corrected  by  simple  combinations  of  linear 
functions  of  distance.  These  massage  functions  were  applied 
to  the  four  forwardmost  angles  and  the  averaged  backward  angles. 
Attempts  were  made  to  keep  angular  errors  in  the  parametrized 
data  compared  to  the  calculations  less  than  20  percent,  however, 
larger  errors  may  occur  in  isolated  cases.  Special  attention  was 
given  to  the  four  forwardmost  angles  since  they  constitute  most 
of  the  importance  of  the  angular  distribution  as  well  as  behave 
in  a manner  that  lends  itself  to  parametrization  by  smooth  func- 
tions. The  averaged  backward  angles  were  not  compared  pointwise 
directly  with  the  data  base  because  at  some  points  the  data  base 
values  contained  zeros  or  exhibited  other  nonphysical  behavior. 

The  averaging  scheme  assures  the  physical  trends  in  the  data. 
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3.2  PROMPT  GAMMA  RAYS 


The  prompt  gamma  ray  data  base  was  generated  in  a manner 
similar  to  that  for  neutrons.  Multigroup  cross  sections  were  pro- 
cessed from  point  data  using  the  SMUG^^^^  code.  DTFXRAY^^^^  cal- 
culations were  performed  for  a source  in  each  of  the  18  gamma 
ray  groups.  The  resulting  energy-angular  fluence  was  then 
processed  to  subtract  the  unscattered  fluence.  The  scattered 
fluence  4)  ^ (E^  , r , E , u)  is  a function  of  source  energy,  E^  ; 
range  r;  detector  energy,  E;  and  angle  cosine  with  respect 
to  the  source  receiver  ray,  u.  The  parametric  modeling  was 
performed  with  the  goal  of  achieving  accuracies  in  modeling 
the  integral  fluence  quantities  such  as  energy  spectra  with- 
out compromising  the  modeling  of  the  more  differential  data. 

The  total  scattered  fluence  function,  F(Eg,r),  was  fit 
with  a function  of  the  form: 


W(p) 

where  p is  the 
the  constants  of 


A,  A. 

= A^  + A2  + A^P  + + —2 

O 

distance  parameter  in  gm/ cm*^ 
the  parametrization . 


and  the  A. 

1 


are 


The  first  step  in  modeling  the  scalar  fluence  was  to 
generate  a "reference  shape"  by  fitting  the  spectra  for  a given 
source  band  at  a range  (several  mfp's  from  the  source)  where  the 
spectra  displayed  the  asymptotic  variation.  Next,  the  energy 
fluences  from  each  mesh  point  of  the  data  base  were  normalized 
to  the  reference  shape  such  that  the  value  of  juxtaposed  fluences 
agreed  at  one  energy  normalization  point,  Ej^. 

Several  different  techniques  were  used  to  generate  the 
spectral  shape  at  a given  distance.  The  scalar  fluence  func- 
tion integrated  over  energy  is  normalized  to  unity  by  defini- 
tion, thus  only  the  relative  variation  of  shape  at  a given 
distance  has  to  be  computed  from  the  fits.  The  reference  shapes 
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for  the  gamma-ray  sources  were  chosen  at  a distance  of  3 mfp  (at 
the  source  energy)  and  modeled  with  a third  degree  polynomial, 

S^(E,p^)  = exp{AQ  + UnE  + A2(ilnE)^  + A^CHnE)^} 

2 

where  is  the  reference  distance  in  gm/cm  . 

A five  parameter  fit  in  distance  was  then  made  of  the  shape 
at  the  detector  energy  which  corresponds  to  the  source  energy,  E^ 

S(Eg,p)  = explBp  + ilnp  + B2(2np)^ 

+ B^(enp)^  + B^(£np)"^} 


The  spectral  shape  for  energies  above  the  normalization  energy 

(0.75  MeV)  was  computed  with  S (E,p  ),  and  S(E  ,p),  by 

o o s 


S(E,p)  = exp 


iln  S^(E,p^) 


!ln  E-£n  E 


N 


\e,n  S(Eg,p)  £n  S^(Eg,p^)j 


E„<E<E 
N — s 


The  shape  below  Ej^  was  assumed  to  be  identical  to  the  re- 
ference shape 


S(E,p)  = S^(E,p^) 


E 


"N 


Then  the  normalized  scalar  fluence  is  given  by 

S’(E,p)  = S(E,p)//  S(E,p)  dE 


Special  techniques  were  used  for  the  lower  source  energy 
groups.  Here,  the  fluence  in  each  group  was  fit  as  a function 
of  distance; 

S^(p)  = exp{AQ  -t-  Aj^  8,np  -t-  A2()lnp)^  A2(*’np) 

2 

p < 142  gm/cm 
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S.(p)  - expftn  S.(U2)  + s.(522)  - tn  S.(142); 

142  < p < 522  gm/cm^ 

where  S^(p)  is  the  unnormalized  shape  for  energy  group  i 
(i  £4),  and  the  quantities  S^(142)  and  S^(522)  were  saved 
in  ATR. 


Comparisons  of  the  parametric  fits  of  the  normalized 
scalar  fluence  and  the  data  base  values  indicated  the  fits 
were  generally  within  107o  of  the  data  base  except  at  the  .511 
MeV  annihilation  peak.  A constant  factor  was  used  to  correct 
this  anomaly. 

The  angle  dependent  ratio  function  is  generated  from  the 
data  base  by  the  equation 


R(Eg,r,E,M)  = 


(Eg,r,E,y) 

J <l>  (E  .r.E.p) 


The  parametrization  of  the  ratio  function  for  photon  transport 
was  considerably  more  difficult  than  fitting  the  scalar  fluence 
function.  The  ratio  function  is  three  dimensional  for  a given 
source  energy,  adding  one  more  independent  variable  to  the  para- 
meterization. Furthermore,  several  of  the  angular  quadrature 
points  had  zero  fluence  values,  particularly  close  to  the  source. 
The  zero  values  arise  from  negative  flux  fix-up  in  DTFXRAY  which 
is  employed  where  the  gradient  of  the  fluence  is  so  steep  that 
extrapolation  to  adjacent  quadrature  points  results  in  a negative 
fluence . 

The  first  techniques  that  were  tried  in  fitting  the  ratio 
function  were  unsuccessful.  These  attempts  avoided  the  zero  fluence 
problems  by  fitting  a cumulative  or  integral  function  defined  by 
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— - — - ' ^ ^ 


I(E^,r,E,p)  -J  ^)(E^,r,E,p)  dy 

The  fluence  could  then  be  retrieved  by  differentiation  of  the 
integral  fit.  This  technique  appeared  at  first  to  provide  a 
smoother  function  to  fit  across  quadrature  points  with  zero 
fluence  values.  Although  the  integral  function  could  be  fit 
with  good  accuracy,  the  functions  differentiated  from  the  fits 
were  not  of  satisfactory  accuracy  when  compared  to  the  data 
base  values. 

The  method  finally  employed  involved  grouping  the 
sixteen  angular  fluence  values  into  a set  of  six  groups.  The 
ratio  function  for  each  of  these  six  angular  groups  was 
m.odeled  with  two  dimensional  fits  in  detector  energy  and 
range  for  each  source  energy.  The  ratio  function  for  the 
si.xteen  angular  bins  was  then  reconstructed  from  the  six 
group  fits. 

The  ratio  functions  for  the  gamma  rays  were  collapsed 
to  the  six  groups  according  to  the  data  presented  in  Table  5- 


TABLE  5.  ANGLE  COSINES  OF  COLLAPSED  ANGULAR 

GROUPS  FOR  GAMMA  RAYS  FROM  AN 

ASYMJ1ETRIC  S,,  SET 
16 


f 

[ Angular  Index 

1 

Angular  Cosine  j 

1 

1 

0.9983100  1 

2 

0.9808660  1 

3 

0.9558195  1 

4 

0.9245530 

5 

0.5174825  ! 

6 

-0.7786290 

The  collapsing  consisted  of  averaging  the  angular  denendent 
fluence  over  the  quadrature  points  in  each  narticular  group. 
The  new  data  xvas  normalized  at  the  lowest  energy  according 
to  the  formulas 
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I 


1 

J 

j 

J 

i 

i 

R.(E.p)  = f^(p)  R^(E,p)  i = 1,6 
r:(e^,p)  = 1 

where  the  source  energy  variable  has  been  dropped  for  con- 
venience, E^  is  the  lowest  detector  energy  in  the  data  base 
(0.02  MeV) , and  the  index  i refers  to  one  of  the  six  collapsed 
angular  groups.  The  functions  R|(E,p)  were  then  fit  with  third 
degree  polynomials  in  range  and  detector  energy  by 

R]^(E,p)  = explApCp)  T A^(p)(inE)  + A,,(p)(tnE)“  + A^(p)(enE)^} 
and 

Aj(p)  = exp{B^Q  + Bj^CHnp)  + Bj.^(Snp)^  + B^^(£np)^}  , j = 0,3 

The  actual  number  of  coefficients  required  was  reduced  by  over 
a factor  of  four  by  combining  coefficients  over  several  groups. 

The  normalization  factors  f^(p)  were  found  to  be  ap- 
proximately constant  in  p.  These  factors  were  simnly  saved 
in  an  18x6  array,  one  for  each  source  and  angular  group. 

3.3  X-RAYS 

The  X-tay  data  base  was  generated  in  the  same  manner  as 
the  prompt  photon  data  base.  The  DTFXRAY  discrete  ordinates 
code  was  used  to  determine  the  X-ray  fluence  for  a source  in 
each  of  the  18  X-rav  groups.  The  parametr izat ion  of  the  X-ray 
transport  data  was  also  performed  in  a manner  similar  to  that 
for  gamma  rays.  However,  the  reference  shapes  for  the  X-ray 

9 

data  base  were  chosen  as  the  shape  at  the  distance  53.8  gm/cm“ 
and  fit  with  a third  order  polynomial : 
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S^(E,p^)  = exp{AQ  + £nE  + A2(ilnE)^  + A^C^-nE)^} 

A fit  of  the  fluence  in  the  fourth  group  ^.25-30  keV)  was  made  as 
a function  of  range 

2 

S(E^,p)  = explB^  + ilnp  + B2(£np)  } 

where  E^  = 27.5  keV  and  for  all  energies  below  Ej^  (50  keV) 
the  spectral  shape  was  computed  by 

S(E,p)  = exp{an  S^(E,p^)  + (In  S(E^,p) 

- £n  S(E^,p^))}  E < Ej^ 

Several  specialized  corrections  were  made  to  shapes  just  below 
Ej^  to  increase  the  accuracy  of  the  models. 

For  energies  above  E^^,  a difference  function  was  formed 
from  the  data  defined  by 

iln  S (E,p  ) - In  S(E,p) 

= en  S^(Eg,p°)  - iln  S(Eg,p) 

This  function  was  fit  both  in  energy  and  range 

D(E,p)  = Aq(p)  + A^(p)  £nE  + A2(p)(£nE)^  + A^CpX^nE)^ 


where 


2 3 

A^(p)  = B^q  + B^^  anp  + B^2(*^^p) 


+ B^^(£np) 


-1 


0 < i < 3 
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Finally,  for  energies  above  the  shape  was  computed  by 


S (E , p)  = exp 


Jin  S (E,p  ) - D(E,p)  • {Jin  S (E  ,p  ) 
o o o s o 


- Jin  S(E^,p) 


E > E. 


where  the  shape  at  the  detector  energy  corresponding  to  the 
source  energy  was  fit  in  range 


S(Eg.p)  = exp{CQ 


+ Cj^p 


+ C2P^  + C 


3 . r 

,p  + 


.p-h 


with  five  parameters. 

Similar  to  the  gamma  ray  models,  special  techniques  were 
used  for  lower  source  groups.  For  the  third  and  fourth  source 
group,  each  non-zero  fluence  group  was  fit  in  range.  The 
second  source  group  was  fit  in  range  with  a fifth  order  poly- 
nomial . 

The  normalized  scalar  fluence  function  for  the  X-ray 
sources  was  calculated  by  normalizing  to  unity  the  spectral 
shape  computed  using  the  techniques  discussed  above. 

The  angular  data  are  treated  similarly  to  the  gamma  rays. 
The  data  are  grouped  with  the  X-ray  data  base  being  collapsed 
into  the  six  groups  according  to  Table  6.  The  four  most  forward 


TABLE  6.  ANGLE  COSINES  OF  COLLAPSED  ANGULAR  GROUPS 
FOR  X-RAYS  FROM  A SYMMETRIC  S^^  SET 


Angular  Index 

Angular  Cosine 

1 

0.98940 

2 

0.94460 

3 

0.86560 

4 

0.75540 

5 

0.36980 

6 

-0.75540 

! 


•i 


\ 
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angles  were  each  allocated  one  group, 
were  normalized  at  the  lowest  energy, 

R.(E.p)  = f.(p)  r;(e,p) 


Again,  the  functions 
Ej^,  by  the  equations 

i = 1,6 


Rj^(Ej^,p)  E 1 


i = 1,6 


and  the  R^(E,p)  functions  were  modeled. 

A 4x6  array  consisting  of  the  R^  for  the  highest  four 
energy  groups  (including  the  source  group)  and  the  six  angle 
groups  was  saved  at  two  ranges.  The  value  of  R'  for  any 
range  for  the  highest  four  energy  groups  and  six  angular  groups 
was  found  by  linear  interpolation  between  the  two  ranges.  For 
any  energy  group  below  the  four  highest  in  energy,  R^  is  found 
by  linear  interpolation  in  energy  between  the  fourth  from  the 
highest  energy  group  and  the  value  at  the  lowest  energy  which 
is  required  to  be  unity  by  the  normalization. 

The  six  angular  fluence  values  are  computed  for  both  the 
gamma  and  X-rays  by  the  following  relationship 


^^(Eg,p,E)  = F(Eg,p)  • S(E^.p,E) 


R|(E,p).f. (p)  , i = 1,6 


A straightforward  interpolation  scheme  as  a function  of  angle 
cosine  is  then  used  to  obtain  the  17  angular  values  of 
4^(E^,p,E).  The  original  data  base  had  a low  energy  cutoff  of 
10  keV , but  for  some  high  altitude  applications  the  lower  energy 
fluence  can  be  important.  The  decision  to  employ  a 10-keV  cut- 
off in  the  earlier  versions  was  based  largely  on  the  results  of 

( 12) 

Banks  and  Coleman.^  Their  calculations  indicate  that  energy 
fluence  from  photons  below  10  keV  drops  at  least  five  orders  of 
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2 

magnitude  in  about  2 gm/cm  of  air--about  10  meters  at  sea 
level  but  2 km  at  40  km  altitude.  The  extension  of  the  photon 
data  base  below  10  keV  was  performed  by  treating  the  uncollided 
component  only.  Table  7 gives  the  group  structure  and  cross 
sections  used  internally  to  determine  the  uncollided  fluence. 
The  18  fine  groups  are  collapsed  to  two  groups  (.1-1.  and  I.-IO 
keV)  for  output. 

TABLE  7.  LOW  ENERGY  X-R.4Y  TOTAL  CROSS  SECTIONS 


Group 


Energy 

Boundaries  (keV) 


Cross  Section 
(cm" 1) 


1 

0.1 

- 

0.2 

24.3608 

2 

0.2 

- 

0.3 

7.29217 

3 

0.3 

- 

0.4 

3.44374 

4 

0.4 

- 

0.5 

6.39955 

5 

0.5 

- 

0.6 

15.5199 

6 

0.6 

- 

0.7 

11.9949 

7 

0.  7 

- 

0.8 

8.35611 

8 

0.8 

- 

0.9 

6.07217 

9 

0.9 

- 

1.0 

4.60307 

10 

1.0 

- 

2.0 

1.54431 

11 

2.0 

- 

3.0 

0. 325121 

12 

3.0 

- 

4.0 

0.  124977 

13 

4.0 

- 

5.0 

6.23234(-2) 

14 

5.0 

- 

6.0 

3.40691(-2) 

15 

6.0 

- 

7.0 

2. 18828(-2) 

16 

7.0 

- 

8.0 

1.42225(-2) 

17 

8.0 

- 

9.0 

9. 54224(-3) 

18 

9.0 

- 

10.0 

6.8529(-3) 

Read  as  6.23239  x 10 
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That  is,  the  18  fine  group  structure  is  used  internally  in 
ATR,  however,  only  two  summary  groups  appear  in  the  output. 
The  two  groups  represent  the  accumulation  of  results  for 
those  energy  limits.  Thus,  if  4>^(r)  represent  the  output 
of  the  response  (fluence  or  dose)  in  the  first  two  groups 
and  stands  for  the  internal  representation  of  the  fine 

group  structure,  then  the  following  relationships  hold: 

-a . r 

^2 

0^(r)  = I 0 (r) 

j=jl 


where 


r 

S. 

J 

V. 

J 

a . 

J 

i 

j 

^1 


the  slant  range 

the  source  spectrum  values 

the  dose  weights  (Table  18)  or  other  fluence  weights 
as  applicable 

the  cross  section  values  listed  in  Table  7 
1 or  2 

an  index  into  the  fine  group  structure 
1 and  ^2  = 9 when  i = 1 

10  and  j2  = 18  when  i = 2 


Section  6.1  deals  with  the  X-ray  source  specification. 
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3.4  FISSION  PRODUCT  GAMMA  RAYS 


The  fission  product  gamma  ray  data  base  consists  of  em- 
pirical equations  which  describe  the  growth  of  the  debris  cloud, 
the  rise  of  the  debris,  the  mass  thickness  of  air  between  the 
debris  and  the  detector,  and  the  transport  of  the  photons 
through  the  disturbed  atmosphere. 

The  air  transport  data  for  the  fission  product  model  was 
generated  by  Version  3 of  ATR.  The  source  values  given  in 
Table  14  are  the  result  of  integrating  the  time  dependent  source 

(13’) 

distribution  of  Engle  and  Fisher^  'over  all  time  and  regrouping 
the  results  into  the  ATR  energy  group  structure. 

The  following  four  dose  responses 

• Henderson  Tissue 

• Concrete  Kerma 

• Air  Kerma 

• Silicon  Kerma 

were  calculated  using  the  dose  response  function  values  of 

Table  16.  These  dose  values  were  calculated  at  several  points 

2 

of  slant  range  as  a function  of  rho-R  (i.e.  units  of  gm/cm  ) 

2 

to  a maximum  distance  of  300  gm/cm  for  an  air  density  of 
-3  3 

1.11  X 10  gm/cm  . The  dose  values  were  then  fit  by  the 
following  six  parameter  equation: 

f(r)  = 4TTR^D^(r)  = exp^a^  + a^  /r  + a2r  + + “4  “3 

where  D. (r)  represents  one  of  the  four  dose  responses  as  a 
^ 2 

function  of  r in  units  of  gm/cm  , R in  cm  and  the  a^  are 
the  parameterization  coefficients.  Thus,  24  coefficients  give 
an  excellent  representation  of  the  transport  data  for  the  four 
dose  responses.  In  fact,  the  relative  error  is  on  the  order 
of  a percent. 


The  source  strength  of  fission  product  gamma  rays  is 
given  by: 

S(T)  = 1.45  X 10^^  FYG(T) 

at  time  T after  fission;  where  it  is  assumed  that  there  are 
23 

1.45  X 10  fissions/KT,  a fission  fraction  F,  a total  yield 
Y in  KT  and  a time  dependent  fission  product  gamma  emission 
rate  of  G(T)  photons/sec/fission.  The  decay  rate  G(T)  is 
given  by 


g 

G(T)  = g'yrp  photons/sec/fission 

which  is  obtained  by  a fit  to  the  Engle  and  Fisher  data.  Note 
that  the  number  of  fissions/KT  is  included  as  a mtiltiplier,  and 
thus  the  number  of  "particles"  generated/KT  must  be  included  as 
a normalization  factor  for  neutrons  and  gamma  rays  if  they  are 
to  be  compared  with  the  fission  product  results. 

Thus  the  total  dose  is: 

D.d.r)  . , 

^ 4ttR‘'(1+.  87T) 

It  should  be  noted  that  r and  R are  related  to  each  other  for  a 
particular  geometry  configuration.  R is  the  line-of- sight  dis- 
tance measured  in  cm  and  r is  the  mass  thickness  of  air  in  that 

2 

same  distance  measured  in  gm/cm  . For  a fixed  R,  r can  change 
due  to  the  shock  wave. 

Further  modifications  such  as  air-ground  interface  effects, 
cloud  rise  and  hydrodynamic  effects  are  applied.  The  air-ground 
interface  correction  factors  utilized  in  the  fission  product  data 
base  are  from  a modified  first-last  collision  model.  Table  8 
presents  a summary  of  air-ground  functions  which  are  used  in 
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TABLE  8. 


AIR-OVER-GROUND  CORRECTION  FUNCTIONS 
FOR  FISSION  PRODUCT  GAMMA  RAYS 


Mean  Free 
Path 


C(Hg)  C(H^) 


0. 

00 

0. 

834 

0. 

875 

0. 

01 

0. 

843 

0. 

882 

0. 

02 

0. 

847 

0. 

887 

0. 

05 

0. 

855 

0. 

897 

0. 

10 

0. 

86  7 

0. 

910 

0. 

25 

0. 

893 

0. 

935 

0. 

50 

0. 

921 

0, 

959 

0. 

75 

0. 

943 

0, 

973 

1. 

00 

0. 

961 

0, 

981 

1. 

50 

0. 

983 

0. 

991 

0 

^ . 

,00 

0. 

992 

0, 

995 

3. 

00 

1. 

000 

1 , 

000 

the  following  manner.  If  H is  the  source  height  and 

the  target  height  (both  measured  relative  to  the  ground),  then 

the  dose  at  slant  range  r*  is  corrected  by  the  following  factors 

Djlr)  = D.(r)  C(H^)  C(H,j.) 

where  C(U,, ) and  C(H.„)  are  obtained  from  the  values  in  Table  8 
S 1 

by  interpo lat  icn . and  H,j,  are  converted  to  a corresponding 

number  of  mean- f ree-paths  by  assuming  an  average  energy  of  1.3 
MeV  for  the  fission  product  gamma  rays. 

The  cloud  rise  effects  are  inherently  time  dependent, 

therefore,  R becomes  t ime  dependent  : R(T)  and  its  formulation 

comes  from  the  fitted  equation  for  the  cloud  height  H„  measured 

(14)  ^ 

from  the  burst  point  as  a function  of  time  T: 

H^(T)=  61  T'^^  meters 


* r is  mass  thickness. 


where  Y is  the  weapon  yield  in  KT.  If  R„  is  the 
horizontal  range  between  the  source  and  detector  , then 
the  following  equation  holds  for  R(T)  : 

R(T)  = + {H^(T)  + Hg  - 

and 

r(T)  = R(T)  p 

— O 

where  p is  the  average  density  in  gm/cm  of  air  between 
and  H^.  The  fireball  rise  model  is  independent  of  height 
of  burst  and  thus  may  differ  from  other  more  detailed  models 
such  as  that  in  the  IDEA^^^^  code. 

A modification  to  the  mass  thickness  r(T)  is  made  to 
account  for  the  hydrodynamic  effect  which  creates  a perturbation 
of  the  air  mass.  The  following  expression  gives  the  scaled  dis- 
tance, J; 


J = 3.473  X 10^  meters 

where  Y is  the  yield  in  KT  as  before , P is  the  air  pressure 
2 

in  dynes/cm  at  the  source  altitude  and  Q is  given  by  the 
following  relationship: 

Q = 1 if  Hg  > 2Rpg 

1 + -O  ('  - If 

where  H„  is  the  source  height  above  the  ground  and  R^^  is 
the  radius  of  the  fireball. 
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Given  the  yield  Y in  KT  and  letting  represent  the  ratio 
of  air  density  at  the  burst  point  to  the  density  at  sea  level, 
the  following  expressions  are  used  to  determine  the  fireball 
radius ; 


R 


A 


.033 


Y 

R 


1/3 

TO 

D 


Rj.  = 20  Und  + .05R^) 


R - ^ 

O ^»Rg  Ra 


R„„  = 1000  r3R  + .1  R (27.5  R-^  - 5 R„)]  (meters) 
FB  L o o o o J 


The  fireball  radius  is  also  used  to  control  a validity  indicator 
on  the  output,  that  is,  an  asterisk  will  appear  on  the  dose  out- 
put corresponding  to  distances  that  are  less  than  two  fireball 
radii  away  from  the  source. 


A scaled  unitless  time  is  calculated  from  the  scaled 
distance  by: 

L = .011831  (P/p)^^^  j 

where  p is  the  air  density  at  the  source  altitude. 

On  the  basis  of  J and  L there  are  two  other  quan- 
tities which  are  calculated: ^ 

If  L < 6.55  then 

C = J .55974  + .951  and 


V = .95134  + .17505 
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If  L > 6.55  then 


C = J(L  + .67279)  and 
V = 1.0  . 

A parameter  X is  obtained  from  V: 

X . - 1) 

7V  + 35 


Finally,  the  modified  mass  thickness. 


r’(T)  is  given  by 


r'  (T) 


r(T) 


21V^  /% 

56V^  - 35  \ ^ 


for  Rg  < C 


r'  (T) 


r(T) 


/ 35(V^-1) 
\56V^  - 35 


for  Rg  > C 


where  Rg  is  the  original  slant  range  between  the  source  and 
detector . 


The  final  equation  for  the  time  dependent  dose  including 
the  cloud  rise  and  hydrodynamic  effect  is  given  by: 


D,(T.r)  - 1.16  X 10^3 
^ 4ttR(T)^(1  + .87T) 


This  must  be  integrated  over  time  in  order  to  obtain  the 
fission  product  gamma  ray  dose; 


D^(r)  = 


1.16 


X 10 

4 7T 


23 


FY 


r 


(T)  1 dT 
R(T)"(1 + .87T) 
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i-s  chosen  to  be  60  seconds  and  the  integration  is  per- 
formed by  a 12-point  gauss  quadrature  in  each  of  three  regions 


1.  0 - 1 seconds 

2.  1-10  seconds 

3.  10-60  seconds 

The  air-ground  interface  correction  factors  are  administered 
for  each  term  of  the  integral  since  the  configuration  of  the 
source  height  is  time  dependent. 
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4.  BUILT-IN  SOURCES,  RESPONSE  FUNCTIONS 
AND  ANGULAR  QUADRATURE 


Although  ATR  has  the  flexibility  to  utilize  any  neutron, 
gamma  ray  or  X-ray  source  or  any  response  function,  it  is  con- 
venient to  have  "standard"  sources  and  responses  that  can  be 
used  with  a minimum  of  input.  All  built-in  sources  and  response 
functions  are  in  the  22  neutron  and  18  gamma  ray  group  structure 
given  in  Tables  9-11.  Note  that  the  neutron  source  energy  range 
is  narrower  than  the  response  energy  range.  It  is  also  noted 
that  the  neutron  group  structure  starts  at  the  lowest  energy  group 
with  the  highest  energy  - 15  MeV  - being  the  upper  boundary  of  the 
last  group.  Similarly,  the  photon  groups  increase  with  energy. 
Table  12  contains  the  group  structure  for  X-rays. 

Two  neutron  source  spectra  are  built-in.  A weaponized 
fission  source  and  a typical  thermonuclear  source  are  con- 
tained in  ATR-4.  The  built-in  values  are  given  in  Table  13. 

The  prompt  and  delayed  gamma  ray  sources  are  given  in  Table  14. 
There  is  no  built-in  X-ray  source  spectrum  although  a black 
body  spectrum  may  be  easily  specified. 

Many  common  response  functions  are  built-in.  Tables  15-17 
give  the  response  functions  for  neutrons,  gamma  rays  and  x-rays. 

As  noted  in  Section  3 , a low  energy  (<10  keV)  X-ray  fluence  is 
calculated  internally.  The  response  functions  used  in  this 
energy  range  are  given  in  Table  18.  For  specific  problems, 
other  response  functions  may  be  utilized,  see  Section  6.4. 

All  sources  and  response  functions  are  isotropic  but  an 
energy-angular  dependent  response  may  be  used  by  performing  the 
integration  with  the  energy-angular  fluence  external  to  the  ATR 
code.  The  angular  quadrature  weight  and  angular  intervals  are 
given  in  Table  19. 
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TABLE  10.  DETECTOR  ENERGY  BOUNDARIES  FOR  NEUTRONS  (MeV) 


TABLE  11.  SOURCE  AND  DETECTOR  ENERGY  BOUNDARIES  FOR 
PROMPT  GAMMA  RAYS,  AND  DETECTOR  ENERGY 
BOUNDARIES  FOR  SECONDARY  GAMMA  RAYS  (MeV) 


1. 

0.02 

- 0.05 

10. 

1.33 

- 1.66 

2. 

0.05 

- 0.1 

11. 

1.66 

- 2.0 

3. 

0.1 

- 0.2 

12. 

2.0 

- 2.5 

4. 

0.2 

- 0.3 

13. 

2.5 

- 3.0 

5. 

0.3 

- 0.4 

14. 

3.0 

- 4.0 

6. 

0.4 

- 0.6 

15. 

4.0 

- 5.0 

7. 

0.6 

- 0.8 

16. 

5.0 

- 6.5 

8. 

0.  8 

- 1.0 

17. 

6.5 

- 8.0 

9. 

1.0 

- 1.33 

18. 

8.0 

- 10.0 

TABLE  12.  SOURCE  AND  DETECTOR  ENERGY  BOUNDARIES  FOR  X-RAYS  (keV) 


TABLE  13.  SOURCE  SPECTRA  FOR  NEUTRONS 


Energy  Boundaries 
(MeV) 

Fission  Weapon 
Source  Values 
(Fraction  in  Group) 

Thermonuclear 
Source  Values 
(Fraction  in  Group) 

1.07(-7)  - 2.9(-5) 

0 

0 

2.9(-5)  - 1.01(-4) 

0 

2.00(-3) 

1.01(-4)  - 5.83(-4) 

0 

2.40(-2) 

5.83(-4)  - 3.35(-3) 

0 

1.22(-1) 

5 

3.35(-3)  - 0.111 

2.227(-l) 

3.65(-l) 

6 

0.111  - 0.55 

1.693(-1) 

1.02(-1) 

7 

0.55  - 1.11 

2. 159(-1) 

8.50(-2) 

8 

1.11  - 1.83 

1.468(-1) 

6.20(-2) 

9 

1.83  - 2.35 

1.060(-1) 

2.80(-2) 

10 

2.35  - 2.46 

5. 743(-3) 

5.00(-3) 

11 

2.46  - 3.01 

2.871(-2) 

1.90(-2) 

12 

3.01  - 4.07 

5.481(-2) 

2.60(-2) 

13 

4.07  - 4.97 

1. 177(-2) 

1. 70(-2) 

14 

4.97  - 6.36 

1. 832(-2) 

1.80(-2) 

15 

6.36  - 8.19 

1.274(-2) 

1.47 (-2) 

16 

8.19  - 10.0 

7.342(-3) 

1.41(-2) 

17 

10.0  - 12.2 

0.0 

2. 56(-2) 

18 

12.2  - 15.0 

0.0 

7.06(-2) 
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TABLE  14.  SOURCE  SPECTRA  FOR  GAJ1MA  RAYS 


roup 

Energy  Boundaries 
(MeV) 

Prompt  Fission 
Source  Values 
(Fraction  in  Group) 

Delayed  Source 
Values  (Fraction 
in  Group) 

1 

0.02 

_ 

0.05 

3.084(-2) 

.0363 

2 

0.05 

- 

0.1 

1. 355(-2) 

.0580 

3 

0.1 

- 

0.2 

8. 164(-2) 

. 1052 

4 

0.2 

- 

0.3 

6.872(-2) 

.0942 

5 

0.3 

- 

0.4 

8.678(-2) 

.0835 

6 

0.4 

- 

0.6 

1.7681(-1) 

.1380 

7 

0.6 

- 

0.8 

1.4017(-1) 

.1090 

8 

0.8 

- 

1.0 

1.0042(-1) 

.0835 

9 

1.0 

- 

1.33  1 

1.0729(-1) 

. 100 

10 

1.33 

- 

1.66 

6. 183(-2) 

.0644 

11 

1.66 

- 

2.0 

3.935(-2) 

.0424 

12 

2.0 

- 

2.5 

3.756(-2) 

.0401 

13 

2.5 

- 

3.0 

2.233(-2) 

.0218 

14 

3.0 

- 

4.0 

2.116(-2) 

.0182 

15 

4.0 

- 

5.0 

7.483(-3) 

.0054 

16 

5.0 

- 

6.5 

3.230(-3) 

.0 

17 

6.5 

- 

8.0 

6.790(-4) 

.0 

18 

8.0 

“ 

10.0 

1. 580(-4) 

.0 
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TABLE  15.  DOSE  RESPONSE  FUNCTIONS  FOR  NEUTRONS  (rad/ (n/cm^) ) 
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6(-9)  A.90(-9) 


TABLE  16.  DOSE  RESPONSE  FUNCTIONS  FOR  GAMMA  RAYS  (rad/ (y/cm^) ) 


Upper 

Energy 

(MeV) 

Henderson 

Tissue 

Concrete 

Kerma 

Air 

Kerma 

Silicon 

Kerma 

1 

h ■■  

0.05 

8. 37(-ll) 

5.90(-10) 

4.40(-ll) 

4.13(-10) 

2 

0.1 

3.90(-ll) 

1. 20(-10) 

9. 75(-ll) 

7.25(-ll) 

3 

0.2 

6.60(-ll) 

8. OO(-ll) 

5.92(-ll) 

7.25(-ll) 

4 

0.3 

1.22(-10) 

1.20(-10) 

1.  II(-IO) 

1. 17(-10) 

5 

0.4 

1. 77(-10) 

1.68(-10) 

1.63(-10) 

1.65(-10) 

6 

0.6 

2. 56(-10) 

2.42(-10) 

2. 38(-10) 

2.37(-10) 

7 

0.8 

3.50(-10) 

3. 30(-10) 

3.26(-10) 

3.28(-10) 

8 

1.0 

4.45(-10) 

4.  lO(-lO) 

4.  lO(-lO) 

4. lO(-lO) 

9 

1.33 

5. 30(-10) 

5.05(-10) 

5.05(-10) 

5.05(-10) 

10 

1.66 

6.44(-10) 

6. 15(-10) 

6. 15(-10) 

6.10(-10) 

11 

2.0 

7. 35(-10) 

7.15(-10) 

7. 13(-10) 

7. 12(-10) 

12 

2.5 

8. 75(-10) 

8.40(-10) 

8. 30(-10) 

8.40(-10) 

13 

3.  0 

1.08(-9) 

9.80(-10) 

9. 52(-10) 

9.85(-10) 

14 

4.0 

1.27(-9) 

1. 18(-9) 

1. 12(-9) 

1.20(-9) 

15 

5.0 

1. 59(-9) 

1.46(-9) 

1.  34(-9) 

1.48(-9) 

16 

6.5 

1. 76(-9) 

1.80(-9) 

1. 60(-9) 

1.83(-9) 

17 

8.0 

2. 07(-9) 

2. 18(-9) 

1.90(-9) 

2.28(-9) 

18 

10.0 

2.42(-9) 

2.65(-9) 

2.24(-9) 

2. 80(-9) 
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TABLE  17.  DOSE  RESPONSE  FUNCTIONS  FOR  X-RAYS 


? 2 
(rad/ (x/cm  ■) ) except  the  gold  dose  which  is  ( (cal/gm) / (x/cm  )) 


Group 

Upper 

Energy 

(KeV)' 

Henderson 

T issue 

Gold  Kernia 

Air  Kerma 

Silicon 

Kerma 

1 

15 

6. 05(-10) 

6. 39(-14) 

4. 68(-10) 

3.  56(-9) 

2 

20 

2. 89(-10) 

7.  32(-14) 

2.04(-10) 

1.73(-9) 

3 

25 

1. 70(-10) 

4.82(-14) 

1. 24(-10) 

1.00(-9) 

4 

30 

1. 12(-10) 

3.44(-14) 

8.09(-ll) 

6.49(-10) 

5 

35 

7.99(-ll) 

2.59(-14) 

5. 78(-ll) 

4.52(-10) 

6 

45 

5.48(-ll) 

1.83(-14) 

4.07(-ll) 

2.94(-10) 

7 

55 

3.97(-ll) 

1.23(-14) 

3.07(-ll) 

1.84(-10) 

8 

65 

3.46(-ll) 

8.96(-15)  1 

2. 77(-ll) 

1.  30(-10) 

9 

75 

3.37(-ll) 

6.80(-15) 

2.  79(-ll) 

l.Ol(-lO) 

10 

90 

3. 54(-ll) 

1.65(-14) 

3.06(-ll) 

8. 19(-11) 

11 

105 

4.06(-ll) 

1.93(-14) 

3. 55(-ll) 

7.  12(-11) 

12 

120 

4. 73(-ll) 

1. 53(-14) 

4. 16(-11) 

6.  78(-ll) 

13 

140 

5. 58(-ll) 

1. 22(-14) 

4.95(-ll) 

6.88(-ll) 

14 

160 

6. 59(-ll) 

9.68(-15) 

5.92(-ll) 

7.  36(-ll) 

15 

190 

7. 96(-ll) 

7.  54(-15) 

7. 19(-11) 

8. 22(-ll) 

16 

220 

9.67(-ll) 

4.83(-15) 

8. 76(-ll) 

9.50(-ll) 

17 

260 

1. 16(-10) 

4.56(-15) 

1.06(-10) 

1.  II(-IO) 

18 

300 

1.40(-10) 

3.61(-15) 

1. 27(-10) 

1.  31(-10) 
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TABLE  i8.  LOW  ENERGY  X-RAY  DOSE  RESPONSE  FUNCTIONS 
(rad/ (x/crn"^) ) except  the  gold  dose  which 
is  ( (cal/g) / (x/cm2) ) 


Group 


Boundaries 

(keV) 


Gold 

Kerma 


Air 

Kerma 

5. 349(-8) 
2.618(-8) 
1. 381(-8) 
1. 339(-7) 
1.279(-7) 
1. 118(-7) 
9.043(-3) 
7.450(-8) 
6. 198(-8) 
3.477(-8) 
1. 189(-8) 
6.083(-9) 
3.935(-9) 
2.601(-9) 
1.898(-9) 
1.418(-9) 

1 . 100(-9) 
8. 773(-10) 


Silicon 

Kerma 


2.083(-7) 
1.  761(-7) 
1. 323(-7) 
1.001(-7) 
6.295(-8) 
4.893(-8) 
3.918(-8) 
3.213(-8) 
2.684(-8) 
2.918(-8) 
7.063(-8) 
3.837(-8) 
2.413(-8) 
1.652(-8) 
1.209(-8) 
9. 310(-9) 
7. 300(-9) 
5. 780(-9) 


TABLE  19.  angular  QUADRATURES,  COSINES  AND  ANGLES 


Index 

Quadrature 

Cosine 

Quadrature 

Weight 

Cosine 

Boundary 

Lower  Upper 

Angle 

Boundary 

Lower  Upper 

1 

-.989401 

1. 35763(-2) 

-1.0 

-.97284 

180.0 

166 . 6 

0 

-.944575 

3. 11268(-2) 

-.97284 

- .91058 

166.6 

155.6 

3 

-.865631 

4. 75793(-2) 

-.91058 

-.81542 

155.6 

144.6 

4 

-.755404 

6.23145(-2) 

-.81542 

-.69080 

144 . 6 

133.7 

5 

- .617876 

7.47980(-2) 

- .69080 

-.54120 

133.7 

122.8  1 

6 

-.458017 

8.45783(-2) 

-.54120 

- . 37200 

122.8 

111.8 

7 

-.281604 

9. 13017(-2) 

-.37200 

- . 18940 

111.8 

100.9  1 

8 

- . 095012 

9. 47253(-2) 

- . 18940 

0.0 

100.9 

90.0 

9 

. 095012 

9.n7253(-2) 

0.0 

.18940 

90.0 

79.1  i 

10 

. 281604 

9. 13017(-2) 

. 18940 

. 37200 

79.1 

68.2  I 

11 

.45801 7 

8.45783(-2) 

. 37200 

. 54120 

68.2 

57.2  ! 

12 

.617876 

7.47980(-2) 

. 54120 

.69080 

57.2 

46.3  1 

13 

. 755404 

6.23145(-2) 

. 69080 

.81542 

46.3 

35.4 

14 

.865631 

4. 75793(-2) 

. 81542 

.91058 

35.4 

24.4 

15 

.944575 

3. 11268(-2) 

.91058 

.97284 

24.4 

13.4 

16 

— -i 

.989401 

1. 35763(-2) 

.97284 

1.0 

13.4 

0.0 

i 
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5.  AIR- GROUND  AND  EXPONENTIAL  AIR  CORRECTION  FACTORS 


The  ATR  data  base  is  a parameterization  of  infinite  air 
results.  Many  problems  of  interest  involve  the  location  of  a 
source  or  detector  near  the  air-ground  interface  or  near  the 
"edge  of  the  atmosphere”.  A rigorous  solution  to  the  treat- 
ment of  these  two  extremes  would  involve  the  generation  of  a 
much  larger  data  base  than  that  required  for  infinite  air. 

Because  such  a data  base  does  not  exist  and  its  generation  was 
beyond  tbe  scope  of  the  present  work,  an  engineering  approach 
was  taken  in  correcting  the  infinite  air  results  to  account 
for  the  air-ground  interface  and  the  exponential  atmosphere. 

The  earlier  versions  of  ATR  utilized  a first-last  collision 
model  which  was  dependent  on  the  source  height  and  detector 
height  in  mean-free-paths . This  model  which  was  developed  for 
neutrons  was  also  used  for  secondary  gamma  rays.  It  was  ex- 
tended for  use  for  prompt  gamma  rays  and  modified  for  use  in 
treating  the  exponential  air  correction.  It  was  realized  that 
the  logical  reasoning  behind  the  appropriateness  of  a first- 
last  collision  model  does  not  hold  for  secondary  gamma  rays 
and  when  a better  data  base  became  available  it  was  incorporated 
into  ATR- 3. 

The  correction  factor  approach  based  on  a limited  data 
base  leads  to  uncertainties  in  the  results  and  the  user  should 
examine  the  applicability  of  ATR  for  each  particular  problem 
based  on  his  accuracy  requirements.  For  all  correction  factors, 
corrections  are  only  made  to  intensity  and  n^  change  in  spectra 
or  angular  distribution  results  from  a correction  for  the  air- 
ground  interface  or  exponential  atmosphere. 


5.1  AIR-OVER-GROUND  FACTORS  FOR  NEUTRONS  AND  SECONDARY  GAMMA 


RAYS 

For  the  neutron  and  secondary  gamma  ray  transport  model 
in  ATR,  air-over-ground  correction  factors  are  based  on  calcu- 
lations utilizing  the  two-dimensional  discrete  ordinates  code 
(16) 

DOT.  ' Several  dose  responses  were  calculated  in  an  air-over- 
ground geometry  for  a series  of  source  altitudes  up  to  300  m. 
The  calculations  were  performed  with  ENDF-IV  cross  sections  in 
the  22  neutron  and  18  photon  group  structure  of  ATR.  The 
air  and  ground  compositions  employed  in  the  calculations  are 
given  in  Table  20  and  only  two  separate  neutron  sources  were 
used  at  each  burst  height  - the  14  MeV  source  and  the  weapon 
fission  source  given  in  Table  21.  The  source  heights  used 
in  the  calculations  were  1,  50,100  and  300  meters. 


TABLE  20.  AIR  AND  GROUND  COMPOSITIONS 
(atoms /b-cm) 


Element 

Air  (p  = 1.22  g/e,) 

3 

Ground  (p  = 1.7  g/cm 

H 

1 

1 

9.7656-3* 

N 

4.0242-5 

0 

1.0697-5 

3.4790-2 

A1 

4.8828-3 

Si 

1.1597-2 

Read  as  9.7656  x 10  ^ 


Correction  factors  for  air-over-ground  radiation  trans- 
port involve  the  ratio  of  the  air-over-ground  result  to  the 
infinite  air  result  at  the  same  slant  range  and  air  density. 
The  tissue  dose  in  infinite  air  was  calculated  by  ATR. 
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TABLE  21.  WEAPON  FISSION  SOURCE 


Group 

Upper  Energy  (MeV) 

Spectrum 

1 

15. 

0.0 

2 

12.2 

0.0 

3 

10.0 

7.342-3* 

4 

8.18 

1.274-2 

5 

6.36 

1.832-2 

6 

4.96 

1.177-2 

7 

4.06 

5.481-2 

8 

3.01 

2.871-2 

9 

2.46 

5.743-3 

10 

2.35 

1.060-1 

11 

1.83 

1.468-1 

12 

1.11 

2.159-1 

13 

0.55 

1.693-1 

14 

0.111 

2.227-1 

15-22 

0.00335 

0.0 

Read  as  7 . 342  x 10 


In  the  air-over-ground  calculations,  for  each  source  (14  MeV  and 
fission)  there  are  four  burst  heights  and  for  each  burst  height, 
four  target  altitudes  were  selected  at  approximately  0.5,  50, 

100  and  400  meters  above  the  ground. 

If  R)  represents  the  air-ground  tissue  dose  as 

0 o T 

a function  of  source  height  (Hg)  , target  height  and  hori- 
zontal range  (R)  , and  D^(Hg  , , R)  represents  the  infinite 

air  ATR  tissue  dose  results  in  a similar  functional  form,  then 
the  correction  factor  C(Hg,Hq,,R)  is 


This  correction  factor  was  parametrized  as  a power  series : 


C / , H , R \ ~ a„  + a,  R + ...  + a r'^  , 

^ Sq  Tq  ) 0 1 n 

That  is,  for  a given  source  height  and  target  height,  the  cor- 
rection factor  is  curve  fitted  as  a function  of  horizontal  range. 
The  number  of  terms  is  4,  5,  or  6 depending  on  which  polynomial 
fit  gave  the  best  approximation.  For  neutrons  and  secondary 
gamma  rays,  two  sources,  four  source  heights  and  four  target 
heights,  a total  of  281  coefficients  were  necessary  to  ade- 
quately represent  the  selected  data  base. 

The  correction  factors  are  applied  differently  depending  on 
the  source  and  detector  heights,  see  Fig.  3, 

Hj,  (100  m)  ■ - ■ » 


Figure  3.  Regions  of  Interest  for  Air-over- 
Ground  Correction  Factors 
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For  a geometry  configuration  in  region  0,  (i.e.  source  less 
than  300  meters  and  detector  less  than  400  meters)  the  cor- 
rection factor  C(Hg,H^,R)  is  evaluated  from  the  parametriza- 
tion.  In  region  5,  C(Hg,H^,R)  = 1.  For  the  other  regions 
the  correction  factor  is  computed  by  the  following  formulae; 

1000-H„ 

1.  = 1 + {C(300,400,R)  - 1}  • ^ 


C(Hg,H^,R)  = 1 + (C^-1) 


1000-H, 


2a.  C(Hg,H.j.,R)  = 1 + { C(300  .H.^,R)  - 1} 


1000-H, 


2b.  C(Hg,H.j,,R)  = 1 + {C(Hg,400,R)  - 1} 


1000-H„ 


200-H. 

3a.  C(Hg.H.j,,R)  = 1 + {C(300  ,H.j,,R)  - 1}  ' 


200-H, 

3b.  C(Hg,H.j.,R)  = 1 + {C(Hg,400,R)  - 1}  • 


1000-H, 

4a.  = 1 + (C(300.H.^,R)  - 1}  • 


C(Hg,H.j,,R)  = + {C(300,H.j.,R)  - C^l  • 


200-H., 


4b.  = 1 + {C(Hg,400,R)  - 1] 


1000-H, 


200-H. 


C(Hg,H.^,R)  = + {C(Hg,400,R)  - Cj^}  • 
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These  calculated  correction  factors  then  multiply  the 
infinite  air  results  in  ATR  to  produce  the  air-ground  corrected 
data.  The  correction  factors  from  the  parametrization  of 
the  14  MeV  source  are  used  for  source  energies  above  5 MeV. 

The  correction  factors  from  the  fission  source  are  used  for 
source  energies  below  5 MeV.  Due  to  the  parametrization  in 
the  different  spatial  regions  some  minor  changes  in  shape  can 
occur  near  the  boundaries  but  these  "perturbations"  are  generally 
less  than  107o.  Since  the  correction  factors  are  based  on  a 
tissue  dose  response,  values  for  other  responses  may  be  in  error. 

5.2  AIR-OVER-GROUND  FACTORS  FOR  PHOTONS 

There  is  no  data  base  or  first-last  collision  model 
similar  to  that  developed  for  neutrons  with  which  to  correct 
for  air-ground  interface  effects  on  the  transport  of  photons 
in  air.  Because  of  the  predominant  forward  scattering  dis- 
tribution of  high  energy  photons,  there  is  less  of  a perturba- 
tion of  infinite  air  results  due  to  the  air-ground  interface. 

The  correction  factors  that  were  employed  for  the  prompt 
gamma-ray  transport  were  based  on  a modification  of  French' s^^^^ 
first-last  collision  method.  Because  the  Klein-Nishina  scatter- 
ing distribution  is  peaked  in  the  forward  directions  for  gamma- 
rays,  those  gamma  rays  which  leave  the  source  in  directions 
pointing  away  from  the  detector  point  (or  target)  will 
have  little  effect  on  the  dose  at  the  detector.  This  is  true 
for  photon  transport  in  both  homogeneous  atmosphere  and  aii - 
over-ground  geometries.  Calculations  bv  Woolson  in  infinite 
homogeneous  air  and  air-over-ground ^ indicate  that  in 

both  cases  the  contribution  to  the  dose  from  backward  directed 
photons  is  generally  about  two  orders  of  magnitude  or  more 
below  the  dose  contribution  from  forward  directed  photons.  This 
suggests  a revision  of  French's  method  to  account  for  the  fact 
that  the  most  important  effect  of  the  ground  is  to  remove  the 
contribution  to  the  dose  from  photons  which  have  left  the  source 
in  directions  toward  the  detector  and  scattered  into  the  ground. 
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The  original  first-last  collision  model  does  not  dis- 
tinguish between  the  directionality  of  the  source  particles 
in  accounting  for  the  effect  of  ground  removal  of  first 
collisions  or  for  the  effect  of  removal  of  the  last  flight 
photons  due  to  the  presence  of  the  ground.  To  account  for 
these  "effects"  the  model  was  modified  for  prompt  gamma  rays 
by  assuming  that  a fraction  of  the  photons  are  unaffected. 

That  is 

F = (1-x)  + xF 

Y 

where  F^  is  the  prompt  gamma-ray  correction  factor,  F is 
the  first-last  collision  model  correction  factor  and  the  para- 
meter X is  a number  between  zero  and  unity.  The  lower  the 
value  of  X the  greater  the  importance  of  directionality,  i.e. 
the  effect  of  ground  removal  of  photons  is  only  important  for 
the  fraction  x of  isotropically  emitted  photons  and  it  is 
assumed  that  (1-x)  photons  are  not  affected  by  the  ground. 

The  value  of  x was  chosen,  rather  arbitrarily,  to  be  .25  and 
does  not  vary  with  photon  energy. 

The  correction  employed  in  ATR  for  the  air-ground  inter- 
face effect  on  prompt  photon  transport  can  be  expressed  by 
the  equation 


())(Hg,H^,R)  = F^(Hg)  F^,(H,j.)  ^ JR) 


where  |)  (Hg , h.^, , R)  is  the  ATR  result  for  source  and  target 
heights  Hg  and  H.^,  respectively  at  range  R and  is 

the  infinite  air  fluence  at  range  R.  The  functions  F^(Hg) 
and  F^(H.p)  are  the  modified  source  and  target  first-last  col- 
lision model  correction  factors  for  prompt  gamma  rays.  These 
factors  are  given  in  Table  22  as  a function  of  the  number  of 
mean- free-paths  from  the  ground. 
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TABLE  22.  AIR-OVER-GROUND  CORRECTION  FACTORS 
FOR  PROMPT  GAMMA  RAYS 


Mean  Free 

Path 

F(H^) 

F(H^) 

0.00 

0.83A 

0.875 

0.01 

0.843 

0.882 

0.02 

0.847 

0.887 

0.05 

0.855 

0.897 

0.10 

0.867 

0.910 

0.25 

0.893 

0.935 

0.50 

0.921 

0.959 

0.75 

0.943 

0.973 

1.00 

0.961 

0.981 

1.50 

0.983 

0.991 

2.00 

0.992 

0.995 

3 . 00 

1.000 

1.000 

This  method  of  correcting  for  the  air-over-ground  per- 
turbations in  gamma  ray  transport  does  not  account  for  the 
increased  dose  over  the  infinite  homogeneous  air  dose  for 
ground  ranges  close  to  the  source  for  source  heights  within, 
say,  100  meters  of  the  ground  which  is  due  to  albedo  scattering 
in  the  ground. 

Correction  factors  for  the  air-over-ground  interface 
effect  for  X-ray  transport  were  derived  from  a series  of  photon 
transport  calculations  performed  by  Coleman  to  generate  EMP 
sources  due  to  the  air-ground  asymmetry.  Coleman  reported 
results  for  1 meter  and  200  meter  source  heights  and  horizontal 
ranges  to  1000  meters. 

An  analysis  of  this  data  indicated  that  doses  at  targets 
above  the  source  point  for  most  cases  were  within  10-207o  of  the 
infinite  medium  doses  for  the  same  slant  range.  This  data  also 
indicated  a general  trend  that  the  dose  at  ground  targets  was 
on  the  order  of  507o  below  the  infinite  medium  dose. 
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A two  dimensional  air-over-ground  correction  function 
C^(R,y)  was  generated  from  Coleman's  data  for  X-ray  transport. 
This  function  is  used  in  the  following  manner 

t>  (Hg.H.^,R)  = C^(R,y)  <J>^(R) 

where  i)  (Hg , , R)  is  the  ATR  result  for  source  and  target 

heights  Hg  and  respectively  at  range  R,  P is  the 

infinite  medium  fluence  at  range  R and  y = H.j,/Hg.  Tabulated 
values  of  the  correction  factor  C^(R,y)  are  given  in  Table  23. 
Functional  values  for  given  R and  y arguments  are  found  by 
interpolation  in  the  table.  Correction  factors  are  employed 
in  this  model  only  if  the  target  height  is  within  200  meters 
of  the  ground. 


TABLE  23.  VALUES  OF  THE  X-RAY  AIR- 
GROUND  CORRECTION  FACTOR 


Range 

(m) 

y H.j,/Hg 

0 

. 5 

1.0 

50 

0.12 

0.201 

.412 

95 

0.16 

0.35 

. 72 

lAO 

0.  17 

0.37 

.85 

185 

0.19 

0.42 

.86 

225 

0.224 

0.45 

. 87 

275 

0.286 

0.54 

.89 

325 

0.39 

0.65 

.92 

375 

0.48 

0.68 

1.0 

425 

0.48 

0.70 

1.0 

475 

0.48 

0.72 

1.0 

525 

0.49 

0.73 

1.0 

575 

0.50 

0.74 

1.0 

645 

0.50 

0.83 

1.0 

715 

0.50 

0.83 

1.0 

785 

0.50 

0.83 

1.0 

855 

0.50 

0.  83 

1.0 

925 

0.55 

0.85 

1.0 

1000 

0.55 

0.90 

1.0 
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5.3 


CORRECTION  FACTORS  FOR  EXPONENTIAL  ATMOSPHERE 


Due  to  the  decrease  in  atmospheric  density  near  the 
"edge"  of  the  atmosphere,  radiation  transport  results  may 
exhibit  a significant  departure  from  results  which  are  scaled 
by  the  mass  thickness  of  air  along  a line  from  the  source  to 
detector.  At  altitudes  above  40  km,  the  dose  may  deviate  sig- 
nificantly from  one  dimensional,  density  scaled  results  because 
of  significant  changes  in  atmospheric  density  over  distances 
comparable  to  a particle  mean-f ree-path . 

Corrections  for  non-uniform  air  density  effects  were  ob- 
tained by  extending  the  first-last  collision  model  to  include 
the  density  variation  versus  altitude  in  the  numerical  integra- 
tion scheme.  The  first  collision  correction  factor  for  source 
altitude  was  obtained  by  numerical  integration  of  the  equation 


FCHj)  . 


/7 


_ -u  p(z) 

W(Z)Xp(Z)e  ^ 

— « j — 

X +z 


+z 


dXdZ 


2 

where  is  the  total  cross  section  in  air  in  cm  /gm,  W(Z)  is 

a weighting  factor  given  by  F(Hg+Z)  to  account  for  multiple 
scattering  and  is  obtained  by  an  iteration  of  the  above  equation. 
A term  to  account  for  reflection  is  zero,  since  neutrons  are 
not  reflected  from  the  top  of  the  atmosphere. 


The  last-collision  correction  for  target  altitude  was  ob- 
tained from  numerical  evaluation  of  the  above  integral  without 
the  weighting  function.  The  numerical  results  for 
F(H.j,)  are  listed  in  Table  24  ns  a function  of  altitude  in  terms 
of  the  mean- free-path  to  the  top  of  the  atmosphere.  These  cor- 
rection factors  were  applied  uniformly  to  the  scaled  infinite- 
air  results  in  ATR  for  both  neutron  and  secondary  gamma  rays. 
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TABLE  24.  EXPONENTIAL  AIR  CORRECTION  FACTOR  FROM 
THE  FIRST-LAST  COLLISION  MODEL  AS  A 
FUNCTION  OF  MFP  FROM  THE  TOP  OF  THE 
ATMOSPHERE 


H(MFP) 

F(Hg) 

F(H.j.) 

0 

0 

0 

0.010 

0.067 

0.180 

0.030 

0.092 

0.229 

0.060 

0.133 

0.300 

0.090 

0.166 

0.  350 

0.140 

0.210 

0.408 

0.210 

0.271 

0.482 

0.320 

0.373 

0.585 

0.490 

0.491 

0.676 

0.770 

0.642 

0.  793 

1.200 

0.808 

0.890 

2.000 

0.920 

0.967 

3.  100 

0.975 

0.980 

4.900 

0.985 

0.990 

7.900 

0.996 

0.998 

H > 7.9 

1.000 

1.000 

Verification  of  the  adequacy  of  this  model  is  not  complete  due 
to  the  lack  of  an  adequate  data  base. 

There  are  no  high  altitude  correction  factors  in  ATR  for 
prompt  gamma  rays  and  for  X-rays.  Results  at  high  altitudes 
for  these  sources  are  infinite  air  results  scaled  to  the  actual 
mass  thickness  of  air.  When  an  adequate  data  base  becomes  avail- 
able, correction  factors  should  be  included  for  photon  sources. 
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In  suininary,  air-ground  correction  factors  based  on 
Henderson  tissue  dose  are  used  to  correct  the  intensity  of 
infinite  air  results  for  neutrons  and  secondary  gamma  rays. 

Only  intensity  corrections  are  made  to  energy-angular  distribu- 
tions and  corrections  are  based  on  results  for  a 14  MeV  source 
(for  source  energies  greater  than  5 MeV)  and  on  results  for  a 
fission  source  (for  source  energies  less  than  5 MeV) . A 
modified  first-last  collision  model  is  used  for  prompt  gamma 
rays  and  a data  base  is  utilized  for  X-rays.  Exponential 
atmosphere  corrections  are  based  on  a modified  first-last 
collision  model  for  neutrons  and  secondary  gamma  rays  and 
no  corrections  are  applied  to  prompt  gamma  ray  and  X-ray  re- 
sults. As  a more  complete  data  base  becomes  available,  both 
the  air-ground  and  exponential  atmosphere  correction  factors 
should  be  reevaluated. 
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ATR  COMT-IANDS 


In  the  development  of  ATR,  attempts  were  made  to  utilize 
a command  structure  which  is  user  oriented,  simple  and  con- 
venient yet  provide  a wide  variety  of  options  for  describing 
the  problem.  The  ATR  command  structure  can  be  divided  into 
four  separate  groups: 

• Source 

• Geometry 

• Output 

• Execut ion 

The  source  commands  deal  with  the  description  of  the  source 
energy  distribution  and  intensity  for  the  various  components 
of  radiation.  The  geometry  commands  are  used  to  describe  the 
burst,  ground  and  detector  locations.  The  output  commands  are 
used  to  insure  that  the  required  quantities  are  available  in 
the  desired  units.  The  execution  commands  provide  the  means  of 
performing  several  calculations  without  repeating  all  the  in- 
put . 

All  of  the  commands  utilize  the  same  basic  structure; 
a typical  command  is  of  the  general  form: 

-''<COMMAND  WORD  - , <UNITS  DEFINITION  ■,  <LIST  OF  VALUES'- 

All  commands  must  begin  with  an  asterisk.  If  an  ATR  command  is 
too  long  to  fit  on  a single  80-character  line,  the  command  may 
be  continued  on  subsequent  lines  not  beginning  with  an  asterisk. 
The  restrictions  on  continuation  are  that  no  single  part  of  an 
ATR  command,  such  as  a number,  may  itself  be  split  into  two 
lines  and  there  should  be  at  least  two  numbers  (values)  on  the 
first  card  image.  The  command  word  is  generally  a mnemonic 
name  for  the  type  of  action  to  be  taken. 
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Figure  4 illustrates  the  various  general  commands  and 
the  portion  of  ATR  affected  by  each.  A list  of  the  command 
words  is  given  in  Table  25  with  a verbal  description  of  the 
command  words  given  in  Table  26,  and  a summary  description 
of  the  ATR  commands  is  given  in  Table  27. 

The  second  field  of  the  general  command,  or  units  de- 
finition, is  delimited  by  commas.  The  units  definition  field, 
which  is  optional,  serves  to  explain  the  meaning  of  the 
numbers  in  the  list  of  values  that  follows.  Typical  units 
include : 


MEV 

KEV 

PER  MEV  - 
COS 


Energy  values 
Energy  values 
Source  values 
Cosine  angle  values 


Specific  sets  of  unit  definitions  are  appropriate  for 
each  command.  A default  unit  definition  has  been  selected 
from  the  set  appropriate  to  each  command,  and  is  used  when- 
ever the  units  definition  with  its  surrounding  commas  are 
omit  ted. 


The  list  of  values  field  of  an  ATR  command  is  used  to 
specify  the  numerical  data  a command  may  require.  A number  in 
a list  may  appear  in  a variety  of  forms  to  suit  the  particular 
user  or  problem.  For  example,  some  of  the  forms  in  which  the 
number  400.  may  appear  are  as  follows: 


400  400.  4.E+2  4. +2  4+2  4000-1 


At  least  one  number  must  appear  in  a list  of  values  element. 
Two  or  more  numbers  are  separated  from  one  another  by  the 
occurrence  of  one  or  more  blank  characters.  Therefore,  the 
user  is  restricted  from  specifying  a number  in  which  internal 
blanks  appear,  or  which  is  split  on  two  or  more  card  images. 

A further  restriction  exists  upon  the  magnitude  of  such 
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TABLE  25.  LIST  OF  COMMANDS 


1.  *z-SOURCE(i)  value 

2.  *z-SVAL,  units,  values 

3.  *z-EVAL,  units,  values 

4.  -'z-NORM  value 

5.  *z-YIELD  value 

6.  *BBODY-y,  units,  values 
7 *FP-y,  value 

8.  *xx,  units,  value (s) 

9.  *Cxx,  units,  value (s) 

10.  *GR0UND,  units,  value 

11.  *EXC,  4PIRSQ 

12.  *DEL  conmiand(s) 

13.  *ST0P 

14.  *F1N 

15.  *TITLE  n 

16.  *PRINT/z/(i^  i2  -•■) 

17.  *WRITE/z/(ij^  i2  ■ • • ) 

18.  *D0SE/z/ 

19.  *CONSTRAINT/z/ (xx  i value),  units 

20.  ''FLUXWT/z/ , units,  values 

21.  *E-REGRP/z/ , units,  values 

22.  *A-REGRP/z/,  units,  values 
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TABLE  26.  SYNOPSIS  OF  ATR  COMMAND  WORDS 


SOURCE (i) 

- 

indicates  that  the  i^  source  option  will 
be  chosen 

SVAL 

- 

indicates  that  an  energy  dependent  source 
intensity  will  be  entered 

EVAL 

- 

indicates  that  the  energy  group  boundaries 
will  follow 

NORM 

- 

source  normalization  value  follows 

YIELD 

- 

yield  of  the  weapon 

BBODY 

- 

indicates  that  a black  body  X-ray  source 
spectrum  will  be  specified 

FP 

- 

indicates  parameters  for  fission  product  source 

XX 

- 

indicates  that  the  geometry  component  and 
associated  values  will  follow 

Cxx 

- 

indicates  the  location  of  constant  response 
calculation 

GROUND 

- 

specifies  the  ground  elevation  above  sea  level 

TITLE 

- 

alphanumeric  character  string  identifying  the 
problem 

PRINT 

- 

specifies  which  response  option  will  be 
printed  as  output 

WRITE 

- 

indicates  which  response  options  will  be 
written  on  an  external  storage  device 

DOSE 

- 

indicates  types  of  dose  to  be  calculated 

CONSTRAINT 

- 

indicates  that  the  coordinates  for  a constant 
response  value  are  to  be  determined 

FLUXWT 

- 

fluence  weighting  factors  which  are  to  be  input 

E-REGRP 

- 

regroups  energy  dependent  fluences  into  desired 
energy  groups 
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TABLE  26.  SYNOPSIS  OF  ATR  COMMAND  WORDS  (Cont'd.) 


A-REGRP 

regroups  angular  dependent  fluences  into 
desired  angular  bins 

EXC 

- 

an  action  command  to  execute  the  problem 
specified  by  prior  input 

DEL 

- 

deletes  individual  geometry  or  output  commands 

STOP 

- 

indicates  the  end  of  a problem,  initializes 
input  flags 

FIN 

_ 

indicates  the  end  of  a session 
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TABLE  27.  ATR  COMMAND  SUMMARY 
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/z/  Denotes  output  particle  type:  Same  as  above  plus  KG  for  secondary  ^aiama  rays,  GG  for 
secondary  plus  primary  p,amma  rays  and  FP  for  dose  values  due  to  fission  product  gamma 
rays  and  is  used  with  DOSE  and  CONSTRAINT.  T indicates  total  results  and  is  used  with 
PRINT,  DOSE  and  CONSTRAINT. 


TABLE  27.  ATR  COMMAND  SUMMARY  (continued) 
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*ST0P  Nullifies  all  previous  input  commands  by  initializing 

input  flags 

*FIN  Signifies  end  of  run.  calls  exit  from  ATR 


numbers  since  a number  is  interpreted  as  a function  of  up  to 
three  integer  parts  (a  whole  part,  a fractional  part,  and  an 
exponent) , none  of  the  parts  of  a number  may  exceed  in  magnitude 
the  greatest  integer  value  appropriate  to  the  host  machine,  nor 
can  the  number  generated  from  these  three  parts  exceed  the  host 
machine's  allowable  range  of  representable  numbers. 

There  are  three  different  numeric  formats  in  which 
data  may  be  input: 

1.  A string  of  numbers  all  separated  by  one  or  more 
blank  characters  from  each  other. 

2.  A starting  and  ending  number  separated  by  an  incre- 
ment specification  in  parentheses. 

3.  A repetition  factor  followed  by  an  asterisk(*) 
character  followed  by  the  number  to  be  repeated. 

Numeric  entries  of  all  three  types  may  be  mixed  and  can  be  used 
in  any  order. 

A detailed  discussion  of  the  input  commands  is  given  in 
the  next  four  sections. 

6.1  SOURCE  COMMANDS 

The  source  commands  are  number  1-7  in  Table  25. 

There  are  several  ways  that  an  ATR  source  spectrum  can  be 
specified.  For  each  particle  in  ATR  there  are  standard, 
built-in  source  spectra  which  can  be  used  or  new  values  may 
be  input.  The  following  specific  commands  may  be  used. 

SOURCE (i)  value 

z = N,  G,  X for  neutrons,  prompt  gamma  rays,  and  X-rays  re- 
spectively, and  has  the  same  meaning  for  other  source 
related  commands. 

i = 1,  2,  3,  4,  .5,  or  6 with  meanings  that  are  described  in 
Table  28  for  each  particle. 
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TABLE  28.  SOURCE  OPTIONS  FOR  ATR 


Value  of  i 
in  *z- SOURCE  (i) 

Neutrons 

Prompt 
Gamma -Rays 

X-Rays 

1 

Fission 

Fission 

Black  Body 

2 

Thermonuclear 

Not  Used 

Not  Used 

3 

12.2-15  MeV 

Not  Used 

Not  Used 

4 

Histogram 

Histogram 

Histogram 

5 

Point 

Point 

Point 

6 

Fission/ 

Fusion  Mix 

Not  Used 

Not  Used 

For  example,  the  command;  "'N- SOURCE  ( 1)  will  select  the 
internal  spectrum  for  the  neutron  fission  source  (see  Table  13 
for  the  actual  values)  which  then  will  stay  in  effect  until  it  is 
replaced  by  another  neutron  source  command.  (Note  that  the 
neutron  detector  energy  spectrum  covers  a wider  range  than 
the  source  spectrum.)  For  X-rays  the  *X-S0URCE(1)  command 
must  be  accompanied  by  *BBODY-y  commands  to  specify  the  para- 
meters of  the  black  body  spectrum. 


Histogram  and  point  values  of  the  source  spectrum  can 
be  entered  by  using  the  ■'''Z-S0URCE(4)  and  the  ''"z-SOURCE  (5) 
commands  respectively  in  conjunction  with  the  ''Z-SVAL  and 
*z-EVAL  commands  which  are  used  to  enter  source  values  and 
corresponding  energy  values.  Only  the  "'z-SVAL  is  needed  if 
the  source  values  being  entered  correspond  to  the  internal 
energy  structure.  The  source  values  must  be  placed  in  se- 
quence of  low  to  high  energy  values.  When  both  the  source 
and  corresponding  energy  values  are  input, then  the  energy 
values  either  correspond  to  histogram  boundary  values  by 
having  one  more  energy  value  than  source  value,  or  they  cor- 
respond to  point  source  values  in  which  case  there  will  be 
the  same  number  of  energy  and  source  values.  Again,  the 
energy  values  must  be  specified  in  ascending  order,  and  the 
source  values  must  correspond  to  the  energy  values. 


n 
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In  the  event  that  external  energy  values  are  specified, 
routines  are  called  by  ATR  to  regroup  the  source  values  to  the 
internal  energy- source  group  structure.  The  user  must  exercise 
caution  when  the  input  spectrum  to  be  regrouped  contains  ex- 
treme gradients  at  the  end  points  or  embedded  zero  source  values. 
The  maximum  number  of  energy  or  source  values  is  arbitrarily 
limited  to  50. 

Source  spectrum  normalization  (in  units  of  particle/KT) 
and  yield  (in  units  of  KT)  can  also  be  affected  by  using  the 
*z-N0RM  and  the  '-’'z-YIELD  commands.  The  spectrum  will  not  be 
normalized  if  *z-N0RM  is  not  specified  (standard  internal 
spectra  are  normalized  to  1) . The  default  yield  is  1 KT.  The 
total  source  output  which  is  printed  with  the  source  spectrum 
configuration  is  the  product  of  the  normalization  and  yield. 

Some  care  must  be  exercised  in  source  specification  when 
running  sequential  problems  with  ATR.  When  both  source  and 
energy  values  are  needed  it  is  not  adequate  in  sequential 
problems  to  specify  only  the  source  values.  Subsequent  problems 
in  which  only  source  values  are  given  will  assume  the  source 
values  apply  to  the  internal  energy  group  structure.  Therefore 
both  the  energy  and  source  values  must  be  specified  every  time 
one  or  the  other  is  to  be  changed. 

It  is  a common  programming  technique  for  a user  to  set 
up  energy  and  source  values  in  one  problem,  changing  only  the 
yield  and/or  the  normalization  for  subsequent  problems.  When 
both  the  yield  and  the  normalization  are  subsequently  specified, 
there  is  no  confusion.  l^en  the  normalization  is  specified 
at  least  once,  and  only  the  yield  is  changed,  ATR  will  always 
normalize  the  source  spectrum  to  the  given  normalization  before 
administering  the  yield.  If  the  yield  is  never  specified,  it 
is  assumed  to  be  unity.  A problem  can  arise  if  the  normalization 
is  never  specified.  Because  ATR  assumes  no  inherent  normalization. 


it  integrates  the  source  and  uses  the  result  as  the  normaliza- 
tion (but  does  not  save  that  number  as  the  normalization  value 
for  subsequent  problems) , and  then  administers  the  yield.  There- 
fore, if  only  the  yield  is  changed  in  the  next  problem,  ATR  will 
integrate  the  source  spectrum  it  used  for  the  previous  problem 
to  arrive  at  a new  normalization.  In  other  words,  the  normaliza- 
tion for  the  next  problem  is  essentially  the  total  output  from 
the  previous  problem.  Therefore  a normalization  should  be 
specified  at  least  once  even  if  the  standard  internal  source 
spectra  are  used. 

The  "value”  is  only  meaningful  for  z=N  and  i=6  in  which 
case  the  value  (between  0 and  1)  represents  the  fraction  of 
fission  neutrons  relative  to  the  total  number  of  neutrons 
(fission  + 12.2  - 15  MeV) . 

Examples : 

1.  *N-S0URCE(2) 

which  selects  the  internal  thermonuclear  S],  .ctrum 
(no  other  commands  are  necessary  for  the  source 
spectrum  specification) . 

2.  *G-S0URCE(1) 

which  selects  the  internal  prompt  gamma  ray  fission 
spectrum. 

3.  *X-S0URCE(4) 

which  means  that  accompanying  *X-SVAL  and  optional 
*X-EVAL  command(s)  specify  an  external  X-ray  histo- 
gram source. 

A.  *N-S0URCE(6)  .5 

which  selects  a neutron  source  of  fission/fusion 
(i.e.,  12.2  - 15  MeV)  mix  each  contributing  one-half 
to  tlie  total  source  spectrum. 
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~''?-SVAL,  units,  values 
z = N,  G or  X. 

units  = PER  MEV,  PER  KEY,  PER  GROUP,  and  defines  the  units  on  the 
source  values.  If  the  units  are  not  specified  then  the 
default  choice  is  set  internally  to  PER  GROUP. 

values=  source  values  corresponding  to  either  the  internal  energy 
structure  or  the  energy  values  read  in  by  the  *z-EVAL 
command.  The  order  is  assumed  to  be  from  low  energy  to 
high  energy.  >. 

For  X-ray  sources  and  when  the  ''X-SVAL  command  is  used  alone 
without  the  -'X-EVAL  command  (see  below)  then  the  source  spectrum 
is  represented  in  a ^ group  structure:  The  two  summary  groups 
followed  by  the  standard  18  groups.  VThen  both  the  *X-SVAL  and 
*X-EVAL  commands  are  used,  the  low  energy  groups  are  regrouped 
into  the  internal  18  fine  groups  and  the  spectrum  above  10  keV 
is  regrouped  into  the  high  energy  18  groups  for  internal  repre- 
sentation. The  output  will  appear  in  the  20  group  structure 
as  discussed  above. 

Examples : 

1.  *N-SVAL  .01  10  33  100  21. 

which  defines  the  neutron  source  values  for  the  low 
five  energy  groups.  Other  groups  are  automatically 
filled  with  zeros. 

2.  *N-SVAL,  PER  GROUP,  .01  10  33  100  20. 

which  has  precisely  the  same  effect  as  Example  1. 

3.  *X-SVAL,  PER  KEV,  .2  .2  1 .5  .33  .1  .1  .1 

which  puts  source  values  into  the  lowest  eight  energy  ^ 

groups  of  the  X-ray  source  spectrum  where  the  first 
two  groups  are  .1-1.  keV  and  1.  - 10  keV. 
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“Z-EVAL,  units,  values 
z = N,  G.  X. 

units  = MEV,  KEV.  The  default  unit  is  MEV  for  neutrons  and 
prompt  gamma  rays  and  KEV  for  X-rays. 

values  = source  energy  boundary  values  (if  option  (4)  was 

chosen)  or  source  energy  point  values  (if  option  (5) 
was  chosen) . The  order  is  assumed  to  be  from  low  to 
high  energy  values. 

Examples : 

1.  *N-EVAL  lE-7  lE-4  2E-2  .11  .56  2 5 12.2  15 

which  defines  the  source  energy  values  in  MeV  for 
neut  rons . 

2.  *X-EVAL,  KEV,  10  50  170  235 

w’hich  defines  the  X-ray  source  energy  values  in  keV 
although  the  unit  designation  is  not  necessary  since 
the  default  unit  definition  is  keV  for  X-rays. 


*z-N0RM  value 


z = N,  G,  X. 


value  = source  normalization  in  particle/ KT.  The  default  value 
of  the  internal  source  spectra  is  1 and  there  is  no  ex- 
plicit unit  definition  for  this  command. 


Examp les : 

1.  -"-N-NORM  2+2  3 

which  normalizes  the  neutron  spectrum  to  the  value  of 
2 -1023  neutrons/KT. 

2.  *G-N0RM  1 

which  normalizes  the  prompt  gamma  ray  source  spectrum 
to  1 gamma  ray/KT. 


I 


I 


I 


I 
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-'-z-YIELD  value 
z = N,  G,  X. 

value  = source  yield  in  KT . The  default  value  of  the  source 
yield  is  1 and  there  is  no  explicit  unit  definition 
for  this  command. 


Examp les : 

1.  *N-YIELD  100 

which  produces  a yield  for  neutrons  of  100  KT . 

2.  -'G-YIELD  .2 

which  produces  a prompt  gamma-ray  yield  of  . 2 KT . 
'^BBODY-y,  units,  value  (s) 


y = one  of  E,  T or  W representing  the  total  energy 

normalization,  black  body  temperature  or  weight 
associated  with  each  black  body  component. 


units  = KEV,  MEV  for  y=T  (default  unit  in  KEV)  and  is  KEV, 

MEV,  KT  (kiloton)  for  y=E  (default  unit  is  KT) . y=W 
indicates  the  fraction  of  energy  in  each  hlack  body 
component . 

value (s)  = total  energy  value  or  temperature  values  or  weight 
values . 


This  command  only  applies  to  X-ray  source  spectra  which 
are  described  by  components  represented  by  a black  body  dis- 
tribution. The  black  body  spectrum  specified  in  terms  of 
frequency,  v,  is  written  in  standard  notation 


U(v) 


Suhv^  1 

hv/KT-1 

C e 
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and  the  total  energy  released  is, 

A v3 


TuCE)  dE  = 8IL(^  f ^ (^) 

O-'  C-^h^  O-'  e^-1  C-^h^ 

where  X = hv/KT. 

If  we  define  the  energy  normalized  spectrum  by: 


g(E)dE  = 

TT^  e^-1 

so  that 

f g(E)  dE  = 1 . 

then  the  correspodning  number  spectrum  is  given  by: 

f(E)dE  = ^ . 

IT  e - 1 

The  total  number  of  photons  emitted  for  an  energy  normalized 
spectrum  is, 


cc 

J 


f(E)  dE  = 


11  ± r jd..  dx 

oJ  e^-l 


It  can  be  shown  that 


oJ  e^-1 


dX  = r,(3)  • !(3)  = 2.40A 


where  t is  the  Riemann  Zeta  function  and  !’  the  gamma  function 
and 


no 

/ f(E)  CIE  - ^ (2.404) 


. 37 
KT 
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so  that  when 


then 


dE  = 1 


f (E)  dE  = . 37/KT 


For  a 

yw^=l,  and 
photons  in 


series  of  black  body  energy  weights  such  that 

for  a total  energy  output  E,j, , then  the  number  of 
energy  group  given  by 


N — 7- 

g 4 T 

^ 7T 


u E ^,/KT. 

W.  ^g+I  1 

I / f’(E)  dE 

i E*'/KT. 

g 1 


E 

where  f'(E)  = E^/e  -1  and  the  total  number  of  photons  is 


W . 


N„  = F_(.37)  I 


KT^ 


Notes  on  black  body  source  spectrum  specification: 

1.  If  no  y=E  command  is  specified , then  the  default 
value  E=1  keV  will  be  used. 

2.  If  only  one  value  is  given  on  the  y=T  command  and 

no  y=W  command  is  present, then  is  set  to  1 

automatically . 

3.  If  more  than  one  value  is  given  by  the  y=T  command, 
then  there  must  be  the  same  number  of  values  given 

by  the  y=W  command. 


Examples : 


1.  --'-X- SOURCE  (1) 

•'BBODY-E  .2 
'■'BBODY-T  2 5 
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y = Y for  defining  the  yield  in  KT 

= F for  defining  the  fission  fraction 


value  = yield  in  KT  or  fission  fraction 

Since  the  fission  product  model  is  source  spectrum  independent, 
these  are  the  only  commands  that  are  needed  . The  fission 
fraction  is  preset  to  one  in  ATR  and  is  the  default  value. 
There  is  no  default  value  for  the  vicld. 


Examp les : 

1.  *FP-Y  100 

*FP - F .5 

represents  a weapon  vield  of  100  KT  and  a fission 
fraction  of  . 5 . 

6 . 2 GEOMETRY 

The  geometry  configuration  for  ATR  is  illustrated  in  Fig. 2 
(pg.  17)  where  the  component  coordinates  are: 

Hg(HS)  = source  altitude 

H.j,(HT)  = target  altitude 
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Rg(RS)  = slant  range 

R^(RH)  = horizontal  range 
0 (AN)  = slant  angle 

Three  consistent  component  specifications  define  a complete 
geometry  configuration  with  respect  to  the  ground  level  (as 
long  as  one  of  them  is  either  HS  or  HT) , and  the  other  two 
components  can  be  calculated  from  the  three.  Therefore,  ATR 
requires  the  specification  of  three  geometry  components,  one 
of  which  may  have  up  to  50  different  values.  The  specified 
output  results  are  then  displayed  for  the  various  geometry 
combinations.  The  format  of  the  input  command  is 

*xx,  units,  values 

where  xx  is  replaced  by  one  of  HS , HT , RS , RH , AN.  The  unit 
options  include  most  of  the  reasonable  units  that  are  ap- 
propriate. Values  may  be  specified  in  a list  separated  by 
blanks  or  in  the  format:  n^  (n)  n^,  signifying  values  rang- 
ing from  nj^  to  n2  in  steps  of  n.  In  this  case  n2  must  be 
arithmetically  greater  than  n^. 

Two  of  the  possible  geometry  configurations  result  in 
ambiguities.  UTien  RH , RS  and  HT  are  specified,  there  is  no 
inherent  information  whether  HS  should  be  placed  above  or 
below  HT.  In  order  to  resolve  the  ambiguity  the  characters 
and  should  be  used  with  the  -'HT  command  to  indicate 

that  HS  should  be  placed  above  or  below  HT  respectively  (e.g. , 
'■HT  - 1000).  The  other  ambiguity  occurs  when  RH , RS  and  HS  are 
specified  and  the  placement  of  the  target  height  is  in  question. 
In  this  case,  the  + and  - characters  are  used  with  the  *HS  com- 
mand to  place  the  target  above  or  below  the  source  respectively. 
If  no  character  is  specified  then  the  default  is  +. 

There  is  also  an  option  in  ATR  to  move  the  ground  to  a 
desired  altitude.  It  is  effected  by  the  *GR0UND  command  and 
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all  specifications  involving  HS  and  HT  are  interpreted  relative 
to  the  ground.  l^ien  sequential  runs  are  computed  with  ATR,  the 
ground  is  not  automatically  reset  to  zero  unless  respecified  by 
another  '-'GROUND  command  or  a '’-'STOP  command  is  encountered. 

To  specify  the  fixed  coordinates  for  a constraint  calcu- 
lation, the  command:  ''-Cxx  should  be  used  in  con-junction  v^7itb 
the  -'CONSTRAINT  command.  In  this  case  xx  has  the  same  meanings 
and  options  as  the  other  geometry  components  (except  RS , the 
slant  range,  is  not  allowed  as  it  is  the  geometric  component 
to  be  determined  by  the  constraint  problem) . 

"XX , units,  value (s) 

XX  = one  of  RH , RS , HT , HS  and  AN  denoting  horizontal 

range,  slant  range,  target  height,  source  height, 
and  slant  angle  respectively  (Fig.  2) . 

units  = one  of  M (meters) , KM,  MILE,  YD,  KFT , and  FT  for 

the  distance  specification  and  one  of  DEG  (degrees) , 
RAD  (radians),  and  COS  (cosine)  for  the  angle  speci- 
fication. Because  of  the  symmetry  of  the  geometry 
configuration  about  the  source  the  range  of  admissible 
angles  is  from  -Q0°  to  -^90°.  Default  unit  for  the 
distance  specification  is  M (meters)  and  for  angle 
it  is  DEG  (degrees) . 

value(s)  = one  or  more  values  of  the  corresponding  geometry 
component.  The  maximum  number  of  values  for  any 
coordinate  is  arbitrarily  fixed  at  50. 

As  mentioned  above,  three  of  the  geometry  commands  define 
a complete  ATR  geometry.  Each  geometry  component  may  have 
several  values  in  principle;  however,  the  intended  use  is  for 
two  of  the  components  to  have  single  values  and  the  third 
component  to  have  one  or  more  running  values.  Another  typical 
use  is  to  have  multiple  values  for  two  components  and  one  value 
for  the  third,  in  which  case  the  two  components  with  the  multi- 
ple values  will  be  successively  paired  and  used  with  the  single 
value  of  the  third  component. 
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There  is  a special  unit  definition  option  (GM)  for 
XX  = RS  which  is  allowed  if  the  following  geometry  configuration 
is  specified;  HS , HT , RS . In  this  case  the  values  associated 
with  the  XX  = RS  command  are  interpreted  as  g/cm  of  the  slant 
range.  The  GM  unit  option  is  restricted  to  this  configuration 
only. 

Examples : 

1.  *HT  +,  KM,  1 

'■-RS.  KFT,  .5  (.5)  10 
*RH  100 

which  specifies  the  target  height  at  1 kilometer*, 
several  values  of  the  slant  range  starting  at  .5 
kilofeet  and  ending  at  10  kilofeet  in  steps  of  .5 
kilofeet,  and  the  horizontal  range  at  100  meters. 

Note  that  the  coordinates  may  be  specified 
in  different  units.  Also,  the  source  altitude  will 
be  calculated  to  be  above  the  target  altitude. 

2.  *HS  2000 

*AN,  DEG,  10  (10)  80 
*RS , KM , 1 

which  specifies  the  source  altitude  at  2000  meters, 
the  slant  angle  ranging  from  10  degrees  to  80  degrees 
in  steps  of  10  degrees,  and  a slant  range  of  1 kilo- 
meter. The  DEG  unit  specification  for  the  slant 
angle  is  superfluous  since  it  is  the  default  unit. 

3.  *H8,  KM,  2 
*HT,  KM,  2 

'■-RS,  GM,  50  120  380  410  550 

which  specifies  coaltitude  source  and  target 

at  2 kilometers  above  the  ground  and  the  slant  range 

of  5 different  values  in  units  of  g/cm2. 

"Cxx,  units,  value (s) 

XX  = RH,  HT,  HS  or  AN. 

units  = same  as  above. 

value(s)  = one  or  more  values  specifying  the  geometry  component. 
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Two  of  these  constraint  geometry  commands  must  accompany 
the  “CONSTRAINT  command  with  one  of  the  geometry  components  hav- 
ing a single  value  and  the  other  having  one  or  more  values.  If 
both  have  more  than  one  value  then  ATR  will  arbitrarily  pick  the 
first  value  of  one  of  the  coordinates  as  the  fixed  value. 


Examples : 

1 . '“'CHS  , 

KM,  10 

■“CAN 

-90  (10) 

90 

which 

defines 

the  fixed  coordinates  for  a constraint 

calculation.  The  source  altitude  is  fixed  at  10  kilo- 
meters and  the  slant  angle  has  several  values  between 
-90  degrees  and  +90  degrees  in  steps  of  10  degrees. 

2.  '“CHT,  KFT,  0 20  50  80  150  130 
"'CAN  0 

which  specifies  six  different  configurations  for  the 
target  altitude  in  kilofeet  and  zero  slant  angle. 

'-'CtROUND,  units,  value 

units  = same  as  distance  units  above. 

value  = distance  value  of  the  ground  relative  to  sea  level  for 
the  ATR  problem  geometry. 


If  this  command  is  not  specified  then  sea  level  will  be 
f used  for  the  ground  level.  If  this  command  is  specified,  cor- 

i responding  HS  or  HT  specifications  are  interpreted  relative  to 

i the  ground  level. 

I Examples : 


1.  ^GROUND,  KM,  .5 


which  specifies  the  ground  to  be  at  half  of  a kilo- 
meter above  sea  level. 

2.  '“GROUND  500 

which  has  the  same  effect  as  example  (1)  since  the 
default  unit  definition  is  meters. 
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6.3 


OUTPUT 


This  section  contains  the  description  of  commands  that 
either  directly  or  indirectly  affect  the  output  of  results  from 
ATR.  By  specifying  several  of  these  options,  many  different 
problems  can  be  run  concurrently. 

There  is  a problem  counter  built  into  ATR  which  is  auto- 
matically incremented  by  one  every  time  an  "EXC  command  is 
encountered;  however,  each  of  these  problems  can  be  assigned 
a title  by  using  the  "'TITLE  command.  The  problem  counter  is 
reset  to  one  following  a *STOP  command. 

The  primary  output  mechanism  is  embodied  in  the  --PRINT 
and  "'-WRITE  commands,  by  which  the  fully  detailed  angle-energy 
output  of  ATR  response  characteristics  can  be  obtained.  It  is 
also  possible  to  limit  the  detail  of  the  output  to  obtain  an 
integral  radiation  response  environment.  Table  29  contains  the 
possible  ATR  response  options  corresponding  to  options  in  the 
*PRINT  and  *WRITE  commands  and  Table  30  further  details  the 
units  associated  with  the  responses.  Of  course,  the  output 
is  labeled  with  the  proper  physical  units  for  all  of  the  responses. 

The  read  and  print  logical  unit  definitions  correspond  to 
the  industry  standards.  The  built-in  logical  unit  values  are; 

5:  input  unit  for  command  inputs 

6:  output  unit  associated  with  "'-PRINT  command 

8;  output  unit  associated  with  "-'-'WRITE  command. 

The  '-'dVRITE  command  is  used  to  place  the  results  from  ATR  on  an 
external  storage  device  such  as  disk  or  magnetic  tape  for  further 
analysis . 

Dose  at  several  geometry  coordinates  can  be  obtained  simply 
by  the  use  of  the  "-'-DOSE  command  and  in  this  case  no  other  output 
command  is  necessary.  There  are  eight  sets  of  built-in  dose  re- 
sponses for  neutrons  and  four  for  photons,  see  Section  ^ for 
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values  of  the  response  functions.  In  the  event  that  the  standard 
dose  responses  are  not  adequate,  the  user  can  specify  responses 
of  his  own  by  the  -'FLUXWT  command.  Up  to  five  sets  of  these  out- 
put commands  can  be  specified:  one  each  for  neutrons,  secondary 
gamma  rays,  prompt  gamma  rays,  combined  secondary  and  prompt 
gamma  rays,  and  X-rays.  These  responses  are  stored  until  deleted 
by  the  *DEL  command  or  *STOP  command  or  the  termination  of  the 
program  by  the  *FIN  command. 

Isodose  contour  values  of  any  of  the  integrated  responses 
can  be  obtained  by  the  use  of  the  "CONSTRAINT  command  which  cal- 
culates the  geometry  configuration  corresponding  to  a predeter- 
mined response  level  for  given  values  of  two  other  geometry 
coordinates . 

The  output  of  differential  results  or  results  integrated 
over  only  one  of  the  independent  variables  can  be  displayed 
according  to  user-desired  energy  or  angle  boundaries  by  using 
the  *E-REGRP  or  “A-REGRP  commands  respectively.  The  output 
will  be  regrouped  from  the  internal  energy  and/or  angle  repre- 
sentation to  the  user  specified  boundaries.  These  two  commands 
can  be  used  individually  or  together.  If  a full  angle-energy 
distribution  is  specified  by  the  output  commands,  then  a two- 
dimensional  regroup  will  take  place  before  the  results  are 
displayed.  The  angle  and/or  energy  spectrum  is  regrouped  over 
the  ranges  specified  by  these  commands. 

It  should  be  emphasized  that  for  any  given  ATR  problem 
any  and  all  of  the  output  commands  can  be  specified  for  anv 
and  all  of  the  particles  provided  that  the  set  of  commands  is 
consistent  and  that  the  proper  source  spectrum  has  been  specified. 
An  ATR  problem  in  this  sense  includes  all  cumulative  specifica- 
tions prior  to  an  *EXC  command.  One  example  of  inconsistent 
output  definition  is  to  specify  *E-REGRP  or  *A-REGRP  without 
specifying  a non-zero  response  index  for  either  a *PRTN'^  or 
‘-''’WRITE  command. 
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6.4  SPECIFIC  OUTPUT  COMMANDS 
'^TITLE  n 

n = up  to  74  characters  used  as  a problem  descriptor  title  to 
identify  the  output.  The  title  will  remain  in  effect  until 
a new  title  is  entered  or  a *STOP  command  is  encountered. 

*PRINT/e/  ^ 


z = N,  G,  X,  NG,  GG,  T representing  neutrons,  prompt 

gamma  rays,  x-rays,  secondary  gamma  rays,  com- 
bined prompt  and  secondary  gamma  rays  and  total 
response  respectively.  The  GG  option  sums  the 
quantities  obtained  for  the  secondary  gamma  rays, 
prompt  gamma  rays  and  fission  product  gamma  rays. 

The  T option  is  used  in  conjunction  with  the  PRINT, 

DOSE  and  CONSTRAINT  commands  to  obtain  the  total 
neutron  plus  total  gamma  ray  results  for  the  fol- 
lowing doses: 

Henderson  Tissue 

Concrete 

Air 

Silicon  (ionizing) 

The  resulting  reponse  value,  of  course,  will  depend 
on  which  source  defintions  are  present.  That  is, 
if  both  neutron  and  prompt  gamma  ray  source  spectra 
are  specified  than  all  three  values  of  neutron, 
secondary  gamma  ray  and  prompt  gamma  ray  results 
are  accumulated.  If  only  the  neutron  source  spec- 
trum is  specified,  then  the  results  for  neutron  and 
secondary  gamma  ray  radiation  are  accumulated.  Table 
31  details  the  possible  input  and  output  configurations. 

ij^,i2,etc.  = integer  values  between  0 and  13  separated  by 
blanks  corresponding  to  the  response  options 
listed  in  Table  29.  Note  that  options  9 through 
12  are  used  only  for  neutrons.  If  a minus  sign 

precedes  the  number,  the  fully  differential 
angle-energy  display  of  the  output  will  be  sup- 
pressed and  only  output  detailing  the  response  as 
a function  of  angle  (integrated  over  energy)  and 
as  a function  of  energy  (integrated  over  angle)  is 
displayed  accompanied  by  the  cumulative  total  re- 
sponse. The  minus  sign  makes  no  sense  preceding 
the  0 option  which  is  the  summary  print  option 


J 
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TABLE  31. 


OUTPUT  POSSIBILITIES  AS  A FUNCTION 
OF  SOURCE  SPECIFICATION 


Dose  or 
Constraint 
Activated  By 

Required 

Source 

Specifications 

Effective 

Output 

N 

N 

N 

G 

G 

G 

X 

X 

X 

FP 

FP 

FP 

NG 

N 

NG 

GG 

N,G 

NG+G 

GG 

N , G , FP 

NG+G+FP 

GG 

N,FP 

NG+FP 

GG 

G,FP 

G+FP 

T 

N 

N+NG 

T 

N.G 

N+NG+G 

T 

N,FP 

N+NG+FP 

T 

N,G,FP 

N+NG+G+FP 

Key: 


N : Neutron 

G:  Prompt  Gamma  Ray 

X:  X-Ray 

FP : Fission  Product  Gamma  Ray 

NG:  Secondary  Gamma  Ray 

GG:  Combined  Gamma  Ray 

T:  Total 
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which  includes  total  integrated  quantities  of 
Henderson  tissue  dose,  ionizing  silicon  dose, 
fluences,  currents,  external  flux  weight  and 
average  energy. 

The  summary  printout  also  includes  forward  and 
backward  quantities  which  are  obtained  by  inte- 
grating over  the  forward  angles  (i.e.,  y > 0)  and 
backward  angles  (i.e.,  p < 0)  respectively.  Any 
number  of  appropriate  integer  options  may  be  in- 
cluded and  their  order  is  immaterial.  Also,  print 
option  commands  are  cumulative  through  successive 
*EXC  commands  unless  turned  off  by  using  the  *DEL 
command.  Fission  product  dose  output  cannot  be 
obtained  using  the  PRINT  cotnmand--the  DOSE  and 
CONSTRAINT  command  must  be  used. 


Examples : 

1.  *PRINT/N/(0  -1  4 -9  13) 

which  specifies  the  output  for  neutrons  of  the 
summary  as  well  as  following  responses:  number 
fluence  (omit  full  energy-angle  detail) , energy 
current  (with  full  energy-angle  detail) , concrete 
dose  (omit  full  detail)  and  external  flux  weight 
(with  full  detail) . 

2.  *PRINT/GG/(0  -2  6) 

which  specifies  the  summary  of  re^.  , number 

current  (omit  full  detail)  and  cone.  ^ dose  re- 
sponse for  combined  secondary  and  prompt  gamma 
rays . 

3.  *PRINT/X/(6) 

which  specifies  the  full  detailed  printout  of  the 
gold  dose  response  for  X-rays. 

4.  *PRINT/T/(0) 

which  specifies  the  summary  of  the  total  responses 
from  all  specified  sources. 


95 


*WRITE/z/(i^  . . .) 


This  command  is  analogous  to  the  *PRINT  command  in  terms 
of  option  interpretation  with  the  exception  that  options  15  and 
16  must  be  specified  if  the  results  from  the  *DOSE  or  *CONSTRAINT 
commands,  respectively,  are  to  be  placed  on  the  external  stor- 
age device.  (Note  that  option  14  is  not  used.)  Options  15  and 
16  are  not  used  with  the  *PRINT  command  because  the  *DOSE  and 
'“CONSTRAINT  commands  automatically  turn  on  the  corresponding 
print  options. 

The  consequences  of  the  *WRITE  options  are  precisely  the 
same  as  those  of  the  *PRINT  options  except  that  they  are 
output  to  a different  output  device. 

Examples : 

1.  WRITE/NG/(2  5 8) 

which  specifies  the  detailed  output  of  number  current, 
Henderson  tissue  dose  and  silicon  dose  for  secondary 
gamma  rays  to  the  external  storage  device. 

2.  WRITE/G/(6  15  16) 

which  is  the  specification  of  detailed  output  of  con- 
crete dose,  the  results  from  the  *DOSE  command  and 
the  results  from  the  ’’"CONSTRAINT  command  for  prompt 
gamma  rays  to  the  external  storage  device. 

’^DOSE/z/ 

z = N,  G,  X,  NG,  GG,  T,  FP 

This  command  automatically  prints  integrated  dose  values 
as  a function  of  the  running  coordinate  specified  by  the  geom- 
etry setup.  There  are  eight  different  dose  responses  for 
neutrons  and  four  different  dose  responses  for  photons. 

The  FP  command  is  used  to  activate  the  fission  product 
model.  The  summary  and  other  print  options  cannot  activate 
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the  fission  product  model  because  there  is  only  dose  informa- 
tion in  the  model  and  also  because  it  does  not  contain  any 
angular  dependence.  Table  31  details  the  possible  output 
configurations  for  the  DOSE  and  CONSTRAINT  commands  as  a 
function  of  source  specification.  This  table  illustrates 
which  sources  must  be  defined  in  order  to  get  the  effective 
output  when  the  "particle"  definition  construction  is  used.  For 
example,  if  the  command  *DOSE/GG/  is  specified  for  the  output 
and  neutron,  prompt  gamma  ray  and  fission  product  sources  are 
specified,  then  the  output  will  be  the  stim  of  secondary  gamma 
ray  plus  prompt  gamma  ray  plus  fission  product  gamma  ray 
radiation. 

Examples : 

1.  *DOSE/N/ 

which  specifies  the  dose  vs.  running  geometry  co- 
ordinate for  neutrons. 

2.  *DOSE/X/ 

which  specifies  the  dose  vs.  running  geometry  co- 
ordinate for  X-rays. 

3.  *DOSE/FP/ 

which  specifies  the  fission  product  dose  at  the 
running  geometry  coordinates. 

4.  *DOSE/T/ 

specifies  that  the  total  dose  from  all  specified 
sources  be  determined  for  the  running  geometry 
coordinate. 

•^'^CONSTRAINT/z/ (xx  i value)  , units 

z = same  as  for  DOSE. 

XX  = RS . 

i = is  the  response  designator  (only  options  1-13  are 
meaningful) . 
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value 


2 

= constant  response  value  of  (non  47rR  ) fluence,  current, 
dose  or  external  flux  weight  for  which  the  constraint 
geometry  is  evaluated.  The  value  cannot  be  a 47rR2  re- 
sponse, since  this  may  be  double  valued. 

units  = same  as  the  standard  distance  lanit  definitions.  This 

unit  definition  designates  the  output  units  of  the  geo- 
metry distance  coordinates.  If  no  units  are  specified, 
the  default  is  meters. 

This  command  combined  with  two  geometry  coordinate  speci- 
fications calculates  the  value  of  coordinate  RS  which  corresponds 
to  the  predetermined  response  level  given  by  the  value.  Of 
course,  care  must  be  taken  to  specify  a consistent  and  valid  set 
of  geometry  configurations. 

Examples : 

1.  *CONSTRAINT/N/(RS  13  3+15),  KFT 

This  combined  with  two  commands  (Cxx)  of  constraint 
geometry  coordinates  specifies  isodose  values  of  the 
slant  range  (RS)  for  neutrons  where  the  response 
value  is  3-10^^  of  the  external  flux  weight.  The 
distance  values  will  be  output  in  units  of  kilofeet. 

2.  *coNSTRAINT/GG/(RS  1 5+6),  MILE 

which  specifies  the  constraint  calculation  of  the 
slant  range  on  the  number  fluence  value  of  5 -lO^ 
for  combined  secondary  and  prompt  gamma  rays.  The 
distance  values  will  be  displayed  in  units  of  miles 
on  the  output . 

3.  *CONSTRAINT/X/(RS  6 2+12) 

which  specifies  the  constraint  calculation  for  x-rays 
on  the  slant  range,  the  constant  response  value  of 
the  gold  dose  is  2 10^2  and  the  distance  coordinates 
will  be  output  in  meters  which  is  the  default  when 
the  units  are  not  specified. 

4.  *CONSTRAINT/FP/(RS  5 lE+5) 

specifies  that  the  distance  RS  in  meters  corresponding 
to  a fission  product  dose  of  10^  rads  is  to  be  deter- 
mined. 
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^‘FLUXWT/z/ , units  , values 


z = N,  G,  X,  NG  (same  as  G) . 

units  = arbitrary,  user  specified  units,  up  to  16  characters 

(not  including  a comma)  that  will  appear  on  the  output. 
If  no  units  are  specified  then  the  commas  should  not 
appear . 

values  = flux  weight  values  to  be  used  for  weighting  the  ATR 
energy  spectrum  and  must  correspond  to  the  internal 
energy  structure  (low  to  high  energy  values)  and  con- 
tain 22  neutron,  18  gamma  ray  and  20  X-ray  numbers 
depending  on  x. 

The  flux  weight  values  will  remain  stored  by  ATR 
until  the  termination  of  the  run  or  until  a *ST0P 
command  is  encountered. 


Examples : 

1.  *FLUXWT/NG/,CM**2  ergs/gram,  5.9-8  1.2-8 

8.-9  1.2-8  1.7-8  2.4-8  3.3-8  4.1-8  5-8  6.15-8 

7.15-8  8.4-8  9.8-8  1.2-7  1.5-7  1.8-7  2.2-7  2.65-7 

this  command  specifies  external  flux  weights  for 
secondary  gamma  rays  corresponding  to  the  internal 
18  detector  energy  groups  in  order  of  low  to  high 
energy  values. 

2.  *FLUXWT/N/,  RADS  SIL.  KERMA,  0 0 

0000000  2-11  3.1-11  3.4-11  4.2-11 

4.7-11  4.8-11  4.8-11  5.1-11  5.5-11  6.6-11  6.4-11 

6.6-11  7.5-11 

which  specifies  external  flux  weights  for  neutrons,  the 
22  values  correspond  to  the  internal  neutron  detector 
energy  structure,  the  9 low  energy  values  being  zero. 


~'^E-REGRP/z/ , units,  values 


z = N,  G,  X,  NG,  GG. 

units  = MEV  , KEV.  The  default  unit  is  MEV  for  neutrons  and 
prompt  gamma  rays  and  KEV  for  X-rays. 
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values  = energy  boundary  values  (low  to  high)  which  will  be 
used  to  regroup  output  results. 


In  versions  1 and  2 or  ATR,  problems  have  been  encountered 
with  the  regroup  function.  Version  3 incorporated  a new  method. 
The  basic  problem  can  be  stated  by  the  following.  Given  a set  of 
{f(x^)}  defined  on  the  set  {x^},  one  needs  a mechanism  to  cal- 
culate the  set  {f(Xj)}  corresponding  to  a given  set  } which 
is  different  from  {x^}.  In  actual  practice  some  of  the  {x^} 

may  be  the  same  as  some  of  the  {x.}. 

> J 

These  data  may  be  in  either  histogram  or  point  form  and 
the  regroup  mechanism  must  account  for  both  types  of  data.  The 
approach  utilized  is  based  on  procedures  that  fit  a third  order 
polynomial  spline  function  based  on  five  adjacent  points.  The 
j consequence  of  the  spline  is  that  the  derivatives  are  also 

matched  at  the  transition  points. 

* The  form  for  the  polynomial  is; 

i 0 3 

' y = Pq  P|^(x-x^)  + p.,(x-xp‘'  + p^Cx-xp 

where 

Pq  = ^1 

Pi  = 

P2  = {3(y2-y2^)/(x2-x^)  - 2t^^  - t2}/(x2-x^) 

2 

P3  " ^^1  ^2  ' 2(y2-yj^)/(x2-x^)  }/(x2-x^) 

In  addition,  the  slope  t at  the  middle  point  is  defined  in 
terms  of  the  slopes  m^  at  the  corresponding  points: 

t = (Im^-m^l  m2  + I m2-m^  I m^)  / ( I m^-m^  I + |m,,-mj^|) 
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This  spline  interpolation  scheme  is  applied  to  the  cumu- 
lative integral  and  the  interpolated  values  are  differentiated. 


Examples : 

1.  *E-REGRP/NG/ ,MEV,  . 1 . 5 . 2 3 5 10 

which  specifies  the  six  (seven  boundary  values) 
groups  to  which  the  output  will  be  regrouped  for 
secondary  gamma  rays.  The  unit  definition  is 
superfluous  since  MeV  is  the  default  unit  for  this 
case . 

2.  *E-REGRP/X/10  45  70  95  150  240  300 

which  specifies  the  energy  boundaries  in  keV  (the 
default  unit)  to  which  the  output  will  be  regrouped 
for  X-rays. 

•’’'A-REGRP/z/ , units,  values 
z = N,  G,  X,  NG,  GG 

units  = DEG  (degrees) , RAD  (radians) , COS  (cosine) . Default 
units  are  degrees. 

values  = angle  boundary  values  (low  to  high  cosine  values) 
which  will  be  used  to  regroup  output  results.  An 
odd  number  of  values  must  be  input. 

Since  the  natural  unit  for  the  angle  dependent  data  in 
ATR  is  the  cosine  of  the  angle,  the  user  must  exercise  some 
care  in  the  treatment  of  this  command.  The  default  unit  is 
degrees  and  the  value  range  is  0-180  degrees.  However,  this 
means  that  if  the  cosine  values  (ranging  from  -1  to  +1)  must 
be  entered  low  to  high  then  the  degree  values  must  be  entered 
high  to  low.  Also  since  degrees  and  radians  are  directly  re- 
lated, if  radians  are  used  for  the  units  then  the  values  must 
occur  high  to  low. 
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Examples : 

1.  *A-REGRP/N/ .COS,  -1  -.5  0 .25  .5  .75  1. 

which  specifies  the  angular  regroup  of  neutron  out- 
put to  the  given  cosine  values. 

2.  *A-REGRP/X/,DEG,  180  120  90  60  45  1C  0 

which  specifies  the  angular  regroup  of  x-ray  output 
in  units  of  degrees  although  the  unit  definition  is 
not  necessary  in  this  case. 

6.5  EXECUTION 

Several  commands  in  ATR  are  concerned  with  the  flow  of 
control  of  the  execution.  Many  ATR  commands  fill  corresponding 
buffer  regions  or  set  appropriate  flags,  and  their  order  is 
immaterial.  However,  commands  described  in  this  section  take 
effect  immediately  after  being  entered. 

When  the  *EXC  command  is  encountered,  ATR  executes  all 
concurrent  problems  and  displays  the  results. 

The  *DEL  command  deletes  selected  input  buffers,  i.e., 
it  clears  the  effects  of  certain  previous  commands  or  group 
of  commands. 

The  *ST0P  command  nullifies  all  previous  commands  and 
clears  the  internal  flags;  it  would  normally  follow  an  -'’EXC. 

-^EXC,  4PIRSQ 

This  command  causes  the  execution  of  the  concurrent  ATR 

problems  and  the  display  of  the  results.  If  the  comma  appears 

after  the  *EXC  then  the  resulting  values  will  be  multiplied 
2 

by  4ttR  where  R is  the  slant  range  in  centimeters. 
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"DEL  conmiand(s) 


conmiand(s)  = ona  or  more  of  the  following  list: 


1. 

GEOM 

2. 

DOSE 

3. 

CONSTRAINT 

4. 

E-REGRP 

5. 

A-REGRP 

6. 

PRINT 

7. 

WRITE 

8. 

FP 

Several  deletions  may  be  effected  at  once  with  a single  *DEL 
specification,  provided  that  the  corresponding  command  words 
are  separated  by  blanks.  There  is  no  command  for  source 
spectrum  deletion  since  the  effects  of  a source  can  be  nullified 
by  deleting  the  corresponding  print  options.  In  any  case,  the 
-'STOP  command  also  resets  all  specifications,  even  those  relating 
to  source,  including  the  ones  that  can  be  selectively  reset  by 
*DEL . 

The  special  command  (GEOM)  can  be  used  to  erase  the  set 
of  buffers  associated  with  the  geometry  coordinates.  When  the 
user  specifies  a set  of  coordinates,  the  components  already 
specified  can  be  changed.  When  it  is  desirable  to  define  the 
geometry  with  a different  set  of  components,  the  *DEL  GEOM 
command  should  be  used  and  the  complete  geometry  specified. 

Examples : 

1.  *DEL  DOSE  CONSTRAINT  FP 

this  command  deletes  the  effects  of  the  dose  and  con- 
straint commands  previously  entered  as  well  as  those 
specifying  the  fission  product  source. 

2.  *DEL  GEOM  A-REGRP  E-REGRP 

this  command  deletes  the  previously  defined  geometry 
configuration  as  well  as  the  angle  and  energy  regroup 
boundary  values. 
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''STOP 


This  command  clears  all  flags  and  internal  buffer  areas 
set  up  by  previous  specifications.  ATR  is  initialized  and 
complete  problem  specification  must  follow  the  STOP  command. 

-’'FIN 

This  command  terminates  the  program  with  a FORTRAN  STOP. 
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7.  ATR  APPLICATION  INFORMATION 


The  ATR  code  has  been  implemented  on  several  different 
computers  and  an  attempt  has  been  made  to  keep  the  code  as 
computer  independent  as  possible.  The  code  is  written  in 
FORTRAN  IV  with  Hollerith  constants  kept  at  a maximum  of  four 
characters  to  the  computer  word.  ATR  versions  have  run  on 
the  following  computers:  Univac  1108,  CDC  6600,  CDC  7600, 

DEC  System  10,  GE  635,  IBM  360,  IBM  370.  The  code  can  be  run 
in  a timesharing  interactive  mode  as  well  as  in  a batch  oriented 
environment . 

Version  4 of  the  ATR  code  takes  about  62K  words  of  core  on 
the  DEC  System  10.  The  execution  time  is  typically  one  second 

or  less  for  most  ATR  problems. 

7 . 1 COMMON  ERRORS 

Most  user  errors  on  input  occur  in  the  specification  of 
the  source.  The  convention  is  that  energy  values  and  corre- 
sponding source  values  are  to  be  entered  from  low  energy  to 
high  energy.  Also,  if  point  values  of  the  source  are  specified, 
then  the  number  of  energy  values  must  match  the  number  of  source 
values.  Similarly,  if  a histogram  source  spectrum  is  specified 
in  a group  structure  other  than  that  used  in  ATR,  then  the 
number  of  energy  values  must  be  one  more  than  the  number  of 
source  values.  If  source  values  are  entered  corresponding  to 
the  internal  energy  structure  then  the  energy  values  need  not 
be  specified  but  only  18  values  are  required  for  either  neutron 
and/or  gamma  ray  sources. 

When  multiple  problems  are  run  and  only  the  yield  com- 
ponent of  the  source  is  changed  from  one  problem  to  the  next, 
then  at  least  once  (preferably  for  the  first  problem)  a 


105 


normalization  should  be  entered  so  that  the  code  will  normalize 
the  source  to  the  normalization  value  for  every  problem  before 
multiplying  by  the  yield. 

A common  geometry  error  occurs  if  an  impossible  configura- 
tion is  specified.  For  example,  the  following  set  of  geometry 
specifications : 

*HS,KM,2 
*HT,KM, 1 
*RS  500 

is  in  error  since  the  slant  range  value  is  too  short  to  reach 
from  the  source  to  the  target. 

The  following  describes  the  error  diagnostics  that  are 
possible  from  ATR  Version  4. 

Error  messages  produced  by  ATR  are  of  two  types; 

1.  errors  in  the  form  of  commands 

2.  errors  in  the  content  of  commands 

The  first  is  associated  with  the  input  phase  and  the  second 
with  the  execution  and  output  phase  of  the  code.  The  next 
two  sections  will  elaborate  on  the  error  messages  and  corre- 
sponding probable  causes.  Error  messages  are  always  preceded 
by  two  asterisks. 

7.2  INPUT  PHASE  ERROR  MESSAGES 

1.  **UNRECOGNIZABLE  CONTROL  COMMAND  < input  characters > 

Which  means  that  the  command  just  entered  is  invalid. 
The  <input  characters>  field  displays  the  first  ten 
command  characters  as  interpreted  by  ATR. 

2.  **MISSPELLED  CONTROL  COMMAND- -INTERPRETED  TO  BE 
<ATR  Command > 

In  certain  instances  ATR  ’’forgives"  the  misspelling 
of  a command  and  displays  the  form  of  the  command  as 
interpreted  by  ATR.  Further  processing  of  the  com- 
mand then  takes  place  beginning  with  the  first  blank 
or  comma  character,  whichever  occurs  first. 


3.  **OPTION  <Option  numbers>  NOT  IMPLEMENTED 

This  signals  that  the  conmiand  was  recognized  but  the 
option  corresponding  to  the  command  is  not  imple- 
mented in  this  version.  This  error  message  should 
not  be  seen  in  Version  4 . 

4.  **ILLEGAL  SOURCE  SPEC 

This  covers  a multitude  of  errors  relating  to  the 
form  of  the  *z-SOURCE(i)  command,  such  as  parenthesis 
left  off  or  in  wrong  place,  or  the  value  of  i does 
not  correspond  to  permissible  values. 

5.  **ILLEGAL  UNITS  SPEC 

Which  reports  errors  encountered  relative  to  the 
”,  units,"  specification  for  source  related  commands. 

6.  **ILLEGAL  NORM  SPEC  --  ASSUME  1.0 

Which  will  appear  if  an  illegal  normalization  command 
is  specified.  This  error  need  not  be  corrected  if 
the  value  of  1.0  is  an  acceptable  source  normalization. 

7.  -’'^ILLEGAL  YIELD  SPEC  --  ASSUME  1.0 

Which  is  very  similar  to  (6)  but  applies  to  source 
yield. 

8.  **PARTICLE  DEFINITION  ERROR 

This  will  occur  if  the  particle  definition  of  an  out- 
put command  is  in  error. 

9.  **UNIT  DEFINITION  ERROR 

Which  reports  errors  encountered  relative  to  the 
", units,”  specification  for  non-source  related 
commands.  Further  processing  of  the  command  begins 
with  the  first  blank  or  comma  character  whichever 
occurs  first. 

10.  **VALUE  PARAMETER  ERROR 

Which  means  that  a value  entered  as  part  of  the 
"value(s)"  field  does  not  correspond  to  acceptable 
value  definition.  This  error  message  does  not  cover 
all  possible  value  errors. 
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11.  **0VERFL0W  OF  NUMBER  OF  VALUES 

This  signals  that  the  number  of  values  entered  by 
the  particular  command  exceeds  the  number  of  values 
allowed  (for  most  commands  this  number  is  50) . 

12.  **VALUE  SEQUENCE  ERROR 

Which  usually  means  that  the  "low-to-high"  condition 
was  violated  in  the  "values"  field;  however,  only 
the  first  two  values  are  examined. 

13.  **ILLEGAL  OUTPUT  SPEC 

Which  means  that  an  error  was  detected  concerning 
the  processing  of  a ''PRINT  or  ''WRITE  command. 

14.  **CONSTRAINT  DEFINITION  ERROR 

Which  means  that  the  form  of  a “'''CONSTRAINT  command 
is  in  error. 

15.  **C0MMAND  DOES  NOT  BELONG  - <command> 

Which  means  that  the  argument  of  a *DEL  command  is 
recognized  as  an  ATR  command  but  is  improper  as  the 
argument  of  the  *DEL  command. 

7 . 3 EXECUTION  AND  OUTPUT  PHASE  ERROR  MESSAGES 


1.  '"''ERROR  IN  <particle>  SOURCE  SPEC 

Which  means  that  the  source  specification  correspond- 
ing to  the  <particle>  is  found  to  be  in  error.  Usually 
this  is  due  to  regroup  errors  in  the  source  specifica- 
tion or  non-corresponding  number  of  source  and  energy 
values  or  the  wrong  choice  of  parameters  in  the 
*z- SOURCE (i)  command. 

2.  **N0  SOURCE  SPEC 

This  error  message  usually  follows  those  of  type  (1) 
above  and  signals  that  all  of  the  source  specifications 
are  found  in  error.  This  message  is  usually  followed 
by.  **EXECUTION  ABANDONED. 
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3.  **GEOMETRY  CONSISTENCY  ERROR 

This  means  that  the  specified  geometry  coordinates 
specified  are  inconsistent,  i.e.,  the  three  coordinates 
chosen  by  the  user  do  not  specify  a consistent  geometry. 

A.  **GEOMETRY  DEFINITION  ERROR 

Which  means  that  the  given  values  of  the  geometry  co- 
ordinates specified  are  in  some  way  in  error,  e.g., 
specifying  a longer  horizontal  range  than  slant  range. 

5.  **CONSTRAINT  FLUX  WEIGHT  DESIGNATOR  ERROR 

Which  means  that  the  response  designator  index  specified 
by  the  ^CONSTRAINT  command  is  in  error. 

6.  **CONSTRAINT  EVALUATION  ERROR  i 

Which  means  that  the  evaluation  of  the  constraint 
problem  for  the  given  constraint  geometry  configuration 
cannot  be  obtained  because  of  one  of  the  following 
reasons  depending  on  the  value  of  i: 

1.  geometry  consistency  error 

2.  geometry  evaluation  error 

3.  number  of  initial  iterations  is  greater  than  10 

4.  number  of  gm/cm2  is  greater  than  1050  during  the 

iteration  process  „„ 

5.  number  of  gm/ cm2  is  less  than  10  during  the 
iteration  process 

6.  no  convergence  after  20  iterations 

7.  **AB0VE  CONSTRAINT  RESULT  MAY  BE  IN  ERROR 

Which  refers  to  the  line  of  the  geometry  configuration 
just  above  this  error  message  and  it  indicates  that 
the  computed  slant  range  is  less  than  the  built  in 
acceptable  value. 

8.  **REGR0UP  ERROR 

Which  signals  that  some  kind  of  error  was  encountered 
when  regrouping  the  input  source  spectrum  or  that  an 
error  was  encountered  on  regrouping  output  results. 

One  possible  source  of  error  is  embedded  source  spect- 
rum values  which  are  zero. 
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9 . **EXECUTION  ABANDONED 

Which  usually  means  that  there  was  not  one  acceptable 
source  spectrum  specification. 

10 . **execution  completed 

This  message  appears  whenever  ATR  has  completed  its 
intended  purpose  and  is  ready  for  more  input  commands. 
It  is  not  an  indication  of  an  error. 
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APPENDIX  A 


FISSION  PRODUCT  GAMMA  RAY  DOSE  COMPARISONS 

The  engineering  model  utilized  in  ATR  to  predict  the  fis- 
sion product  dose  is  based  on  empirical  equations  which  provide 
a reasonable  fit  to  experimental  data  from  large  yield  weapons. 
The  parametrized  models  used  in  other  codes  such  as  IDEA^^^^ 
are  generally  much  more  complex  and  are  based  on  detailed 
phenomenology  calculations.  To  provide  a comparison  between 
the  two  models  several  calculations  have  been  made. 

Figure  A-1  shows  the  fireball  rise  rate  for  a 10  KT  burst 
as  it  is  modeled  in  ATR  (independent  of  height  of  burst)  com- 
pared with  the  burst  height  dependent  data  for  IDEA.  Since 
the  phenomenology  used  in  IDEA  (based  on  the  LAMB^^^^  code  from 
the  Air  Force  Weapons  Laboratory)  is  still  undergoing  develop- 
ment, the  more  simplified  model  is  used  in  ATR.  As  further 
data  on  fireball  rise  becomes  available,  the  model  used  in  ATR 
should  be  reevaluated. 

Comparisons  of  dose  values  generated  by  the  fission 
product  model  in  ATR  are  made  with  other  available  data. 

Figures  A-2  to  A-4  present  comparisons  with  curves  generated 
by  the  IDEA  code.  For  these  figures  the  ground  level  in  ATR 
was  placed  at  300  meters  in  order  to  correspond  to  the  air 
density  at  which  the  IDEA  curves  were  generated.  The  target 
height  is  1.75  meters.  Figures  A-5  and  A-6  also  have  the  same 
ground  and  target  parameters  and  they  illustrate  the  source 
and  fireball  rise  model  sensitivity  respectively. 

Figures  A- 7 through  A- 10  represent  other  examples  of  com- 
parisons with  IDEA  and  also  TDM  which  is  a code  used  to  calculate 

(22) 

tissue  dose  in  air-over-ground  geometries.'  For  these  com- 
parisons the  ground  is  at  zero  as  is  the  target  height. 
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In  general  these  comparisons  show  that  the  ATR  model 
represents  the  fission  product  component  of  radiation  reasonably 
well.  The  areas  of  serious  disagreement  with  other  codes  are 
due  to  the  simple  fireball  rise  model  used  in  ATR  since  the 
fission  product  model  is  relatively  sensitive  to  the  fireball 
rise  rate.  The  objective  of  implementing  a fission  product 
model  in  ATR  that  is  very  fast,  with  each  problem  taking  only  a 
fraction  of  a second,  has  been  met.  The  limitations  must  be 
realized  however. 


IIA 
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Figure  A-2.  Fission  Product  Gamma  Ray  Dose  from  a 5 KT 
Burst  100  m Above  Ground 
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Figure  A- 3.  Fission  Product  Gamma  Ray  Dose  from  a 200  KT 
Burst  350  m Above  Ground 


117 


DOSE  ( rads) 


118 


Dose  (rads) 


1000 


1500 


ector 


(in) 


duct 
for  ; 


Henderson  Tissue  Dose 
1 KT  Burst  at  60  Meters 


L22 


1500 


Figure  A- 


TOlKr 

Burst  Detector  Range  (m) 


Fission  Product  Henderson  Tissue  Dose  Com- 
parison for  a 10  KT  Burst  at  5 Meters 


L23 


A 


Burst  Detector  Range  (m) 


Figure  A-10.  Fission  Product  Henderson  Tissue  Dose  Com 
parison  for  a 100  KT  Burst  at  500  Meters 


APPENDIX  B 


SAMPLE  PROBLEMS 

The  following  sample  problems  are  a representative  sam- 
ple of  some  of  the  capabilities  of  ATR  Version  4.  Each  problem 
was  chosen  to  point  out  some  of  the  major  features  of  the  code. 
Each  problem  will  be  described  in  some  detail,  and  the  input 
as  well  as  the  output  is  presented.  If  these  problems  are  run 
in  a single  ATR  session,  then  the  *STOP  command  should  be  used 
after  each  problem  definition. 
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Sample  Problem  1 


This  problem  presents  the  use  of  the  fission  product  part 
of  ATR.  The  fission  product  source  has  a yield  of  10  KT  and 
a fission  fraction  of  0.5.  The  source  and  target  are  both  at 
1 km  height  and  the  separation  slant  range  starts  at  500  meters 
in  increments  of  100  meters  to  a maximum  of  1000  meters.  The 
dose  is  evaluated  for  each  of  the  geometry  configurations 
and  presented. 
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ATR  SAMPLE  PROBLEM  1 
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ATR  SAMPLE  PROBLEM  1 
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Sample  Problem  2 


This  problem  demonstrates  the  total  dose  option  of  the  code. 
The  source  spectra  used  for  the  calculation  are  specified  as 
the  internal  standard  option  of  fission  source  for  both 
neutrons  and  prompt  gamma  rays . Th  normalization  and 
weapon  yield  values  are  not  specified,  thus  the  internal  de- 
fault values  are  used.  The  geometry  configuration  is  such  that 
the  source  is  at  a height  of  1 km,  the  detector  is  at  a height 
of  2 m and  the  horizontal  range  varies  from  100  m to  500  m in 
steps  of  100  m.  For  each  one  of  the  geometry  configurations 
the  total  dose  is  evaluated  and  displayed.  The  total  dose  in 
this  case  is  the  sum  of  the  dose  contributions  from  neutrons, 

secondary  gamma  rays  and  prompt  gamma  rays.  The  dose  values  are 
2 

4itR  quantities. 

NOTE:  If  fission  product  gamma  rays  were  also  included,  they 

would  dominate  the  dose,  since  the  source  consisted  of 
one  neutron  and  one  gamma  ray  but  the  fission  product 
results  include  the  number  of  fissions/KT. 


129 


1 


ATR  SAMPLE  PROBLEM  2 
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ATR  SAMPLE  PROBLEM  2 
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ATR  SAMPLE  PROBLEM  2 - OUTPUT  (Cont'd.) 
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Sample  Problem  3 


This  problem  illustrates  the  X-ray  related  capabilities 
of  the  code.  The  source  definition  is  a composite  of  two 
black  body  spectra,  one  with  a temperature  of  10  keV  and  the 
other  of  20  keV.  The  weighting  of  each  of  the  black  body 
spectra  is  the  same:  0.5.  The  total  output  is  0.4  KT.  The 
geometry  configuration  consists  of  both  the  source  and  the 
target  being  on  the  ground  and  two  slant  ranges  are  used,  one 
at  100  m and  the  other  at  200  m.  For  each  of  the  geometry 
configurations,  the  radiation  summary  is  displayed  on  the  output. 
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ATR  SAMPLE  PROBLEM  3 - INPUT 
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ATR  SAMPLE  PROBLEM  3 
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ATR  SAMPLE  PROBLEM  3 - OUTPUT  (Cont'd.) 
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This  sample  problem  illustrates  the  use  of  the  ^CONSTRAINT 

command.  The  prompt  gamma  ray  source  spectrum  uses  is  the  stan- 

2 A- 

dard  fission  source  with  a normalization  value  of  10  gammas/KT 

and  a vield  of  20  KT . The  neutron  source  is  the  built-in 

23 

thermonuclear  spectrum  with  a normalization  of  6 x 10  neu- 
trons/KT  and  a weapon  yield  of  100  KT.  One  fixed  component 
of  the  geometry  is  the  source  height  at  1 km  the  other  is  the 
slant  angle  which  varies  from  -10  degrees  to  30  degrees  in 
steps  of  10  degrees.  For  each  of  the  five  geometry  configu- 
rations the  slant  range  is  evaluated  for  which  the  response 
value  is  the  same  as  specified  by  the  -'CONSTRAINT  commands.  The 
prompt  gamma  ray  total  fluence  is  constrained  to  have  a value 
of  10^,  the  neutron  total  fluence  is  constrained  to  have  a 
value  of  lO'^  and  the  secondary  gamma  ray  Henderson  tissue  dose 
constraint  value  is  0.1  rads. 
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ATR  SAMPLE  PROBLEM  4 - INPUT 
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ATR  SAMPLE  PROBLEM  4 
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AIR  P«nP|FM  M m'IFR  1 SAmPiF  I’ROHi  Fm  R - ATR  VTRSi  P I. 

»»»<•*»♦*♦»*»**«»*****•*«*«»*•«**»♦ 


fJFilT  Rlir. 

ri'NSIRATill  F(|P 

Si  Al.  I HAIAl,r  AA|(1 

amoiHFP  FliiFaiCF 

r 1 . l OI  +0-1 

AHiTHT/,  HAVGF 

SLAS't  HASGF 

A T AHGF  T All. 

S'liiS'l.F  alt. 

SI  AM  a^'CLF 

S.?3e'>t  + (»3 

S.3104F.  AOi  FI 

7.h290F  +0  1 

1 .OOOOF:  + 03 

-10.000  OFG 

5.S775F+03 

S.S77SfcFn3  M 

1 . OOOOf  +0  3 

T ,OOiiOF  + n3 

o.uon  nFG 

S.7O0«F+03 

S.i'2V3F  *03  f 

2 . 0 1 2 3F  + 0 3 

1 .00(1uF+O3 

1 0.000  i)Ct; 

S.7312F  +03 

b.o>7<>i)t  + n3  M 

3.nn()nF  +03 

1 .OOOoF  + 0^ 

20.1100  OFG 

5.SSl‘TF+03 

6.<iin8E  + o3  II 

O.^Ob'JF  +03 

1 .OOOOF  + 03 

30.000  DtG 

• CALLUlATFI)  FRflM  llTMtR  COlTROlNATFS 
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ATR  SAMPLE  PROBLEM  4 


OUTPUT  (Cont’d.) 


*******  * 

****** 

* * 

A 

******** 

*******( 

k * 

* 

4TR  PKOl3L^‘'  1 

SAmPlF  PWOl'LtM 

a - Aio  s(tPSlji> 

i U. 

******** 

****** 

* * 

* 

******** 

*******( 

k * 

* 

RA  *«A 

c ji'S rwA T>  I 

riip 

SI 

A-,T  4AWf;F  A'Jll 

f'llr-HF  p F 1 HI  Mf  F 

= 1 

, 1 lit  fllS 

♦ HOP  I/!.  WA^i^;^ 

«;LA  l|  r■>A•l(,^ 

* r APC.F,  T Al,  1 . 

SOilPCF  Al  1 . 

SL 

A'  1 A'  S 

LF. 

S«So68f-'  +0^ 

G.»'^P^F-fOi  ‘ 

1 . 8 /4  3 1 F ♦ 0 1 

1 . Olio  Of  ♦ 0 3 

- 

1 1,  . (II  . .1 

ore. 

S.BPH^FtO  ^ 

S.WMH.U+oS  f 

1 . 11  DO  (It  *0  ' 

1 . ODD  01  + 0 S 

il  . 0 0 D 

OF  c; 

6.0800F*Oi 

h.  I S5bt  *T_A  ( 

?.()bHS[  >03 

1 . 0 0 0 1.  f ♦ 0 3 

1 0 . 11  .0(1 

OFG 

h.Ohytf  ♦!>  ^ 

3.P0  /ttt03 

1 , odOiiF  ♦<>  3 

£>'l  , I’Od 

OF 

s.e'iHPF  + ns 

8."t  ■A('K*o.S  ■• 

0 , 0 OS  3(  + 0 3 

1 , n II 0 1 f < 0 s 

S'l  . lllKi 

1 1 F (- 

*r  A|  CHl  A TM.  iM 

Hf 

c. 

iUV!)  I NA  Tt  S 

******** 

****** 

* * 

* 

******** 

»**«**♦ 

k A 

* 

4 T R P K n 1 f-  1 \j ' M ‘-i »-  >0  1 

S A S|P(  t Dl'ilO)  r 

■■1  _ A I P V 1 P S 1 1 1 

'J  . 

******** 

****** 

A * 

A 

******** 

**♦»*»* 

k h 

♦ 

NfinWIlN  GA*mA 

CM^.S  TKA  U I 

»■  1 •U 

St 

AijI  PArgRF  AmO 

HF  'Il'F  PSIU.  1 1 SS|. 

■ = I 

.1  OF-01 

Ahiiwl^,  PA‘jr;r 

S 1 i [ >V  A J 1 . h 

*TAWr,Fi  Al  T. 

SililF'CF  Al  I . 

Si 

i . ( A 1.1 

l.t 

*0  A 

'i.si  r'n  ♦'is  • 

P.lS-i7f  ♦OP 

1 .00, lilt  ♦(!  3 

- 

1 I. . i;'i(i 

OF  1. 

« .6  /O'U  +0'^ 

u . 7 .1  -tF  ♦ (1  i (■ 

\ .0  00  OF  ♦(!  3 

1 . oiMu  F ♦ 0 3 

Il  . II  0 (. 

OF  ,i 

«.76(i2F  ♦0'^ 

■1  . S^bF  ♦O  S ♦ 

1 .H3«3F  + 0 3 

t . o,in(iF  *0  3 

1 11  . ( I'O 

Off. 

U.  7''*S<F  + 0 ^ 

c;.  0(,  7 3F  ♦' -^ 

^.7  1 ph>-  +03 

1 . 0 0 0 0 1 ♦ o ^ 

p,..fi  0(1 

DFi; 

u.o<)^;f  + o^ 

^ , 1 b 7Hl  ♦ fj  S 

3. SO  SOF  f i S 

I . ooOi.F  *0  S 

3 1 . 0 •!  (1 

OF.  (S 

*r4(  c ii 

A I ‘ 1 ■ 1 y ( , .1  , 1 1 K 

******** 

****** 

* * 

* 

A*  AAA-AA* 

*♦••*** 

« A 

* 

♦ ♦EXtC'T  ^ ff  I' 
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This  problem  illustrates  the  output  of  spectral  information 

for  neutrons  and  secondary  gamma  rays.  The  neutron  source  is 

the  built-in  14  Mev  spectrum  with  a normalization  value  of 
7 1 

2 X 10  neutrons/KT  and  the  yield  is  10  KT . The  source  is 

located  on  the  ground,  the  target  is  at  1 km  and  the  slant 

2 

range  is  200  gm/cm  . The  output  details  the  summary  and  spectral 
fluence  values  distributed  over  angle  and  over  energy  but  not 
the  full  angle-energy  distribution. 


i 


ATR  SAMPLE  PROBLEM  5 


* T I 1 1 ^ r-  I ‘ , - A I 

* f ^ 1 1 1 I I I 

-r  I - I ' ,>  s 

-Y  f A I.  . 1 

*M.S  il 
*t*T  ,t"’,  ) 

+ P S » U I'*  » ' 

*p;  1 N W"  Y I ■'  - 1 1 

* i-’  -'  !*'■  I / '■'  I-.  / I - I I 

*f  ’(r  f 


- INPUT 


/f  V ^ 1 T'  i 
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ATR  SAMPLE  PROBLEM  5 - OUTPUT 


ATK  PROBLLH  Mli'UiKk  1 SAMPlK  S - AlH  VEPSld’i  (I. 

1ti(**1t***»»1i******H**-At*A*********** 

NFDTRON  SiliiRCF  T^T^P^4|  VtiFV 

4Jt)R^AL  I Za  1 1 n‘i=?.  iKKit  f ■f'lip'i'j  r rti  i/  = i .i'Ri'f  fOi  -r 

TdTAL  IUnPUI  = Z.nnOF*P  I --H'llk.u' 


S >MPf  r SPF  c TRIM 


ENFRuVfMtV  ) 

N 

N/Mt  V 

FNe  t<l,T  ( ■'‘E  V 1 

1/  '( 

v' 

1 .f/F-n'j-Z.'Ylf 

-r.s 

i' . iinF  - ii  1 

<■  .iinE-i)  1 

? . 3Sf  0 0 - r . 'I  Ell  -fi  Ml 

0 . iiyf  -11 ) 

",  . 1 I-  t 

-■  t 

? . 9f|F  -oS-  1 . II  I ^ 

. 1 . A ( f - y 1 

1 . (!  (lE  ”1)  t 

P . '1  SE  ♦ iMl  - ( . ,1  1 E 

0.  OiiE  -y ) 

y.  .;,1 

J 

1 . 0 1 ^ -ny-B. fi'^f 

“C.l 

•1 . H'F-'.  1 

( .IlOE  -0  1 

A . (1  1 E + iH)  - u . 0 El  4 1'  11 

y . iiyE  -II 1 

y . yt 

- 1 

S.BjF-oy-i. ^s^ 

-0^ 

1 . OliT’-P  1 

0 . y E - 0 1 

'J  . 0 '/ 1 ♦ p .1-  y . R yl  E II  y 

y . y yE  - (1 1 

) . 1,  yi 

- C.  1 

^.3bF-05-1 . 1 IF 

-.M 

0 , <1 1)  T - A 1 

0 . n y E - 0 1 

y . R 7 E 4 y 0 - R . .A  El  E + 11 ,1 

y . Oil!  -11 1 

y . r y E 

-1, 

1 . 1 1 K-0 l-5.Sof 

-f*  ) 

0 . ni)*-  -11 1 

0 .not -0 1 

El . AfiE  + y (,1  -f<  . 1 ”r  + (1  y 

il  . IiyE  - y 1 

y . I 'll 

-V  X 

S.Sc*’  "0  1 - 1 . 1 U + "11 

n . '1  (i  E - 0 1 

1 . i1  0 (■  - 0 1 

«.  ('JEnEIO-I  .'lyE'  + iii 

II  . y (1 K - ii  1 

y . ■ 1 E 

1 

1 . I lF  + OP-1 .R  U 

+ :i  n 

11 . iliiT  -ii  1 

1'  .Ill'll  -11 ) 

1 . iiyf  1)  1 - 1 .R,’l  ♦ I1 1 

y . yyE -y  1 

.1  .1  OE 

1 

1 .H3F+00-«'.3^t 

f'tO 

0 . n(iT  - n 1 

" . noE  -I)  1 

1 .PPE  401-1  .‘Mil  4i'I 

P . y 1 1 F 4 P y 

7.1  1 1 

if  ^ 

*«***«« 

♦ * 

* * * * * 

* * * * 

*4**II4II*1 

* * * * 

* * * 

* 

ATR  F’HnHLtu  MlivirtFw 

1 

SAmPi  F PR'IMl'  '1  5 - 

A 1 P vT  PS  II  '1.  i . 

it  *****  n 

* * 

* * * * * 

* * * * 

»*«**»*«* 

< * * * 

* * -k 

n 

’ GROIINO  ItVfl. 

D.ono'"'', 

■1 . 'I0i)i 

o.nooyFl 

•M  1 1-  S 

AHORIZ.  RA'iFf 

Rr'r 

1 

1 7 Ij  . '1 11  1 1' 

■^ /(.'>**  7> , y . S fi  y K E 1 

, 0 , iWxl  Ml  1 1 h S 

SLAfJl  Range 

RS  = 

1 . / 1 3^’  , 

Pun . TOO! 

s.M'^'eei 

, I e S 

TARl.tT  alt. 

Ml  _ 

1 

' 1 El  . T ‘M  1 

I/Ci**?,  A.P.hIEF  1 

. y.'iPi 

•inis 

sniiRCE  Alt. 

RSr 

■I.OOl’K',, 

'1 . " 1,1  (1  ( 

■i/r  :■* » ?,  y . y n'vr  i 

, ,1  . llOHM  1 1 1 S 

♦slant  angle 

A ■.= 

is.  /I  p.  1 

uRfFS  (LI 

S=  il  . E 1 1 <»  / ) 

*r  ai.c'Jla  » Fo  FRi.<  iiiiEP  r.'i  IPO  I na  tf  s 


********* 

* ♦ * 

* 

* • 

* * * 

* * * * * 

* * * * 

* 

« * 

* * 

Nf  IIT  WIIEJ 

z Y 

1 1 

■!(.tn  1 . 

T'lr  At. 

il  J J 1 IliiM 

KA(  E 

K ■ ) 

1 K"  .V  A >v  P 

tissue  nnsE  (pah 

S*Cii**? 

) 

t . 

M Al  + 1 11 

S.  1 /FI  1 A 

. y Ri.  411  A 

1 

.S/E 

4 E A 

^ . 

»>i)t-  *■  1 ^ 

SlLICilN  DIISE  {RAM 

S*C“**P 

) 

i. 

(1  .AE  E y 0 

1 .SSF 1 1 R 

s 

. 11  ( ♦ ,iR 

? 

. ■|^F 

♦ t 1 

1 . 

AOE  + 1 R 

nueihfk  1 1 uf.'.'ri  ( 

•1 

) 

A. 

'lAE  + 1 M 

y.'liiFiRR 

1 

.y  At  4iMI 

. 7"' 

4RR 

^ . 

/ 1E  4RR 

NUMMFK  (llRRf'ilC 

i 

) 

i. 

'l9E  ♦ 1 N 

y . I'.AE  E P 1 

1 

. 9yf  4 11  A- 

1 

. AAi 

I RR 

1 1 E ♦ R R 

FnFEGY  ElurMCEn 

■i-  V 

1 

y . 

y (1E  4 1 y 

y . yuL  ♦ p.i 

h 

.R  AE  EyR 

K 

..M<E 

' Rl 

A . 

1 lE  ERR 

CNF  WE’ T 

V-Mf-  V 

E 

9. 

/RE  + 1 ■> 

1 . SSF  »RR 

s 

. 'M)f  4 iR- 

1 

.PP' 

4RJ 

f • 

R/E  4RR 

Ext.  e 1 X . M T . ( 

1 

0. 

yyi -y  i 

0 . (lOF  - 0 1 

(E 

.0  01 -y  1 

0 

.yyE 

-y  1 

A ^ 

y Of  -n  1 

average  ENFPtnC’FV 

1 

1 . 

AEiE  * n 1 

y .>«F-oi 

-■1 

.yRE -OR 

• PSE 

-0  1 

Slf-iM 

*««****«* 

* A A 

4 

* * 

* * * 

* * * * * 

4 * » * 

« 

* » 

t • 
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ATR  SAt'lPLE  PROBLEM  5 


OUTPUT  (Cont'd) 


AT«  PHOHLt'l  ■'JIJMHCW  1 SAMi'l.f  S - AT><  VfUSIiiN  'I. 

***«**•*«*****«*******«***•«**«**•* 

I' I S I I 'IT  1 1 in  ,)vF  W A^'G|  t 


ANC, 

OF  GPF  F S 

r> 

'SI 

41- 

S|I-1  '1VtW  1 

'it  5 1 

; Y 

(.H'  1 0 A 1 1 ii'E 

( ' 

'/SU  47 

1 

( 

9 

1 

1 8 0 . n 

- 

1 h 

- 1 . 0 0 0 P 

- 

-II  . 9 /?b 

5.9S0F 

421 

b . 7 594  421/ 

2 

1 8b  , 6 

- 

iss-^ 

-P.97?b 

- 

- P . 'M  0 8 

3.9804 

421 

2.2524  4 21 

i 

1 5b  . b 

- 

1 fJ  ■(  . h 

- b . 9 1 0 r 

- 

-0.«15n 

9 . 0 OSF 

421 

•4. 8b  14  4 ?1 

n 

1 ‘1  ■'!  . b 

- 

13  3.7 

-P . 8 1 S J 

- 

- 0 , b 9 b -1 

9.  1 3 34 

421 

7 . 4 >'  7 4 4 2 1 

s 

13  3.7 

- 

122.'' 

- 1 . 8 9 0 b 

- 

- 0 . S .1  1 ? 

9 . 2 4 8' 

421 

1 . t ' *22 

6 

122.8 

- 

1 n .8 

-0.5412 

- 

-0 . 372' 

9 . 3 41 5 4 

421 

1 . bSol  ♦ 22 

7 

111.8 

- 

101'.  9 

-0.  3 7?(; 

- 

- (•  . 1 8 9 .1 

9.5074 

421 

2. 1 7r4 ♦ 22 

8 

100.9 

- 

90 . 0 

-n . 1 «9-. 

- 

n . 11(1 0 0 

9 . 7 304 

421 

2.  7 394  ^ 22 

9i).n 

- 

79.1 

P . 0 0 0 i; 

“ 

0 . 1 b 9 u 

9.9  301 

421 

3.3284  422 

10 

79.1 

- 

bb.,’ 

<1  . 1 h 9 SI 

- 

0.3  72.1 

5. 1 024 

421 

3.“1 bt  < 22 

1 1 

bp.? 

- 

57.? 

0 . 3 7 2(1 

- 

0.5012 

5.95  34 

421 

o . /j9bl  4 22 

1? 

57.2 

- 

9t>.  3 

).S  i12 

- 

1 1 . r>  9 n '-I 

5 .8954 

421 

b.0U94  422 

13 

«b.3 

- 

55.9 

,1  . b 9 p b 

- 

4>  . 8 1 bo 

b.  3 7 74 

♦21 

b.5ot4  4 22 

1« 

35.^4 

- 

29 . 9 

0 . 8 1 S 9 

- 

0.9108 

7.0184 

421 

5.9884  4 22 

15 

29. u 

- 

15.9 

(1 . 9 1 bo 

- 

u . 9 / 2 0 

7.8n8t  421 

b . 2 7 ii  4 4 22 

16 

13.9 

• 

9 . 0 

0.'-l7?b 

1 .Op  ,1  o 

9.0  1 34  421 

b . <1  2 7 ' 4 2 r’ 

* * 

* * 

* A 

* * * 

* « * » * 

A 

* 41  * * 

A A A A A A A 

A A 

A 

* * 4 4 * 4 

n I S I K J In.I  IM,,  iWfW  hNfhJGV 


EAiFF  GY  ( 'IF  V 

1 S' 

''  ('1/49  A'U.L4 

r'Oiiil  A 1 1 V4  ?"* 

1 /'U  V 

1 

( 1 

1 

0.004-01-9 

. 1 j r - p 7 

9. 7211  42  7 

1 .9bS4  4 2 1 

2 

9. 19t-o  7-1 

. 1 2'  -.14 

3.0554  42  7 

9.1  1 2F  * 2| 

3 

1 . 12f--0b-  3 

. ObF  - ,18 

1 . 3294  42  7 

b.b9l  4 42  1 

9 

3. 0 84 -Ob- 1 

. 9 7 4 - b s 

9 . 7 1 9 4 4 2 b 

1 .0294  422 

S 

1 . 0 74-95-2 

. 9nr -.15 

1 .bSl 4 42b 

1 . 53  1 4 422 

6 

2.9o4-  -05-1 

. 9 1 F - 1 II 

5 . 3„94  42b 

1 . 7 1 81 ( 22 

7 

1 .lllF-O'J-5 

. 9 3r-bo 

1 . 1 824  4 2b 

2.2  /84  4 22 

8 

5.8  3F -09- 3 

. 35r-..  3 

2. 1 1 94  42  1 

2.8b  (4  4 22 

9 

3. 354 -9  3-1 

. 1 1 ' -'M 

1.3314423 

9.29  /4  422 

1 0 

1 . 1 tfc-0 1 -5 

.S„‘--oi 

1.9371  4 22 

5.19  /F  422 

1 1 

S.SoF-01-1 

.111401 

8.9294  421 

5.8  1 94  4 22 

12 

1 . 11  F400- I 

,8  3 4 4 0.' 

3.  35b4  4 2 1 

S . 8b9F  4 22 

13 

1 .8  34  4 0 0-2 

. 354  4 O'. 

2.2914  421 

5, 978442 2 

19 

2.3SF  401-2 

. 9 bji  4 ,1 

2.2('9l  42  1 

8.00 11 +22 

15 

2.984  400-3 

. 9 1 P 4 Op 

1 .2724  42  1 

b . 0 7 1 4 4 22 

1 6 

3.niF  40  0-9 

, 0 74  4 -1 M 

8 . 3 794,  429 

6 . 1 58F  4-22 

1 7 

9.0  /4  400-9 

.9/F  4 09 

8.59  34  4 20 

h.21  34  4 22 

18 

. 9 .9  7f  490-b 

. 38F  *n  1 

b.2924  42(1 

b.  3(10' 4 32 

19 

b.36fc400-8 

. 1 94  4(1  .) 

2.9884  420 

b.  3S9F  4 22 

20 

8, 19E400-1 

. OOF  ♦ ')  1 

1 .buSF  420 

b.  3824  4 22 

21 

l.OOF+Ol-l 

.22r+oi 

9.0b8E4l9 

b.a02F4?2 

22 

1 .22F401-1 

.50F40 1 

1 .31  (14  420 

6.9  394  4 22 

.S  u '• 
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ATR  SAMPLE  PROBLEM  5 


OUTPUT  (Cont'd.) 


******* 

****** 

A A 

* 

* * * * * 

* * * * * 

* * * * 

* * * * * 

Sl'NNARY 

unrnLL , 

Tin  Ai 

HU  Il.lMiP 

l"  AC  OR  AR'l 

F URWARO 

TISSUE  OdSF 

(.»  M'S*  €"»<•? 

) 

i> , ni>F  -0  1 

T. 76r * 1 T 

0 

. nnt -0 1 

? , St^l  * 1 ? 

i.F  Tf  + 1 ^ 

sii  icoM  niisr 

{t-  ^('S*C.^■•  **? 

) 

it.iTfiF-iM 

U.1?E-Hi 

0 

. 0(iF.  - 0 1 

i.  7 /F  •»  1 r’ 

■i.  7UF.  + 1 ^ 

N(.lrtHE«  flUfljO 

( r. 

) 

II  , ,iOF  -n  1 

0 

. C'.F  -11  1 

R.  TiRF  * PP 

/.cor  +;>? 

NDMhEW  rUWKr'iI 

( i; 

) 

0 . (Kll  -'M 

*?? 

0 

.cm  - 01- 

1 . >M  4 ;v 

EKEF(.Y  FI  iIFnCF 

( G-  "F  V 

1 

U . 'I'lt  “‘1  I 

I . nSF  ♦F’lT 

0 

. Clfl  -0  I 

'I.HPI  4PI 

1 . 11 1 F 

FFiFKGY  CUWPF'JI 

( (.-•<(  V 

) 

0 . iif’F  - n 1 

«.R1 F 

0 

, C III  -0  1- 

?.  1 U * ,->1 

R. I PL*?? 

EXT.  E L X . -M  . 

f 

) 

n . lint  -0  1 

1' . (.uF-fi  1 

1, 

. C III  -0  1 

0 . (1  m - T 1 

0 . <1  Of  -d  1 

AVF«f.F-E  FNFKGY 

( f F V ) 

0 . I)I'F  -11  1 

1 . n /F  * no 

0 

.cm  -111 

1 . c 71  - 0 1 

1 . .1  'JF  * 1)11 

***•«<* 

****** 

A « 

* 

* * * * * 

« A A A A 

* « * * 

* * * * * 

ATR  PRIlBLf  M MUMUFR  1 

SAl'RlF 

Ph'MUl.FM  S 

- ATR  VI 

R F 1 1 1'j  II  . 

It  * -k  * n * -k 

A A 

* 

***** 

A A A A A 

A A A A 

* * * * * 

ff 

'lEt'l  F 

('  I S I K I r-i.'  I 1 1 !'■  I H'f  t-  ^ 


ANG 

PF  GREFS 

C 

f'y 

S'."'  IIVFK  Lil?cY 

f i':-ci  1 I ; 1 

( ( 

./STfP  T ( 

1, 

1 

1 Ho  . 0 

- 
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